NASA Technical Memorandum 83135 

NASA-TM-83135 19810024423 

Users’ Manual for a Computer Program 
To Calculate Discrete Frequency Noise 
of Conventional and Advanced Propellers 

R. M. Martin and F. Farassat 
AUGUST 1981 





NASA Technical Memorandum 83135 


Users’ Manual for a Computer Program 
To Calculate Discrete Frequency Noise 
of Conventional and Advanced Propellers 


R. M. Martin and F. Farassat 
Langley Research Center 
Hampton, Virginia 


rVI/NSA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Branch 


1981 



CONTENTS 


SUMMARY 1 

INTRODUCTION 1 

SYMBOLS AND ABBREVIATIONS 2 

THEORY AND EQUATIONS 4 

Theoretical Formulation 4 

Numerical Approach 6 

PROGRAM DESCRIPTION 7 

Main Program - PROPFAN 8 

Main Subroutines 8 

Subsonic Logic 9 

Supersonic Logic 9 

Interpolation Routines 10 

OPERATING INSTRUCTIONS 10 

Input Functions 10 

Namelist Parameters 12 

Output Identification 15 

Output Data 16 

CONCLUDING REMARKS 16 

REFERENCES . . 17 

FIGURES 18 

APPENDIX A - CALCULATION OF EMISSION TIME 21 

APPENDIX B - EXAMPLE CONTROL CARD DECK 25 

APPENDIX C - EXAMPLES OF PRINTED OUTPUT 29 

APPENDIX D - EXAMPLES OF LANGLEY GRAPHICS OUTPUT 55 


iii 



SUMMARY 


This report serves as a users* manual for a computer program for calcula- 
tion of discrete frequency noise of conventional and advanced propellers* The 
computer program was developed at the Langley Research Center by F. Farassat and 
Paul A. Nystrom, The structure of the program and the subroutines describing 
the input functions are discussed in this manual. A detailed definition of 
input variables and their default values and the variables in the output data 
sheet are included here. 

Two versions of the program are available. These differ only in the 
graphic output capability. One version has only printed output capability. A 
second version with extensive graphic output capability is available for the 
computer system at Langley. 

This manual includes four detailed examples of both the printed and graphic 
outputs. These examples may be reproduced by users to check their code on their 
computer system. 


INTRODUCTION 

During the last 30 years, the application of the propeller as a propulsion 
system for large commercial and military aircraft has been neglected. Flight 
speed limitations, excessive noise, and high operating costs are a few of the 
reasons for this neglect. Over the past decade, however, the energy crisis has 
created a new incentive for the aircraft industry to reconsider propellers for 
future aircraft designs- It is well-known that the efficiency of propellers is 
generally higher than that of jet-driven propulsors, but the problem of the high 
noise levels generated by conventional high-speed propellers is a major obstacle 
to their increased utilization. 

An attempt has been made in recent years to develop an advanced blade geom- 
etry to decrease propeller operating noise and to optimize performance. These 
new high-speed blades are designed with substantial out-of-plane sweep and blade 
twist. Model tests have demonstrated the effectiveness of these new designs in 
reducing noise and maintaining high efficiency. 

Since model and full-scale tests of propellers are very expensive and time 
consuming, considerable effort has been made to predict the noise theoretically 
to complement the aerodynamic design calculations. A computer program called 
PROPFAN was developed at the Langley Research Center. This program has the 
capability of predicting the noise generated by a propeller with advanced blade 
geometry operating at both subsonic and supersonic tip speeds. Special effort 
was made during the development of PROPFAN to describe the blade shape as accu- 
rately as possible and to minimize mathematical approximations. The program 
input format also allows the user a great deal of flexibility and convenience 
when applying the code to different engineering problems. 



The first section of this paper describes briefly the theoretical formula- 
tion of this noise prediction method. The method is basically a numerical 
finite-difference approach using two forms of the solution of the Ffowcs 
Williams-Hawkings wave equation. The computing algorithm and a discussion of 
the details of some of the methodology are also included in this section. 

The computer program PROPFAN itself is discussed in the second section. 

The program operating environment, structure, and organization are explained. 
The final section of this paper discusses the operating instructions for exe- 
cuting PROPFAN. The input parameters and functions are defined, and the output 
data are explained. 

Appendix A contains an explanation of the two methods used for calculating 
the emission time. Appendix B contains an example of control cards for a 
multiple-case job as used at the Langley Research Center. Several examples of 
the input functions and printed output are given in appendix C. Examples of 
the graphics output of the version of the program currently operational at 
Langley are presented in appendix D. 
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SYMBOLS AND ABBREVIATIONS 


blade passage frequency 

speed of sound in undisturbed medium 

distance of observer from X 3 ~axis (see fig. Al) 

radial distance along 

= E2/R 

= 0 equation of blade surface 


- t + — = 0 equation of collapsing sphere centered at observer 

c , , 

position X 

indices of summation, i and j = 1, 2, or 3 
function defined by equation (A5) 
leading edge of blade 

local force (per unit area) on fluid in direction i 

helical Mach number 

Mach number in radiation direction 

overall sound pressure level, dB (re 20 yPa) 
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PCA 


pitch change axis 


p’ (x,t) 
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SPL 


TE 
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Po 


acoustic pressure 

nondimensional distance from leading edge in chordwise direction 
based on local chord 

blade outer radius 

length of radiation vector, |x - y| 

-1 • • 

radiation vector, x - y 

unit radiation vector 

surface area of blade 

sound pressure level, dB (re 20 yPa) 

trailing edge of blade 

observer time 

local normal velocity of blade surface 
local velocity of blade surface 

observer position in frame fixed with respect to undisturbed medium 
or in moving frame fixed to aircraft 

observer position in fixed frame 

observer position in moving frame 

source position 

azimuthal coordinate of observer 

curve of intersection of blade surface and collapsing sphere 
Dirac delta function 

rotating frame coordinates; T| is distance from origin in polar 
coordinates (see fig . A1 ) 

angle between blade surface normal and radiation vector 
density of undisturbed medium 
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T 


source time 


T* emission time 

(f) azimuthal coordinate of source at time T > 0 

ip azimuthal coordinate of source at T = 0 

0) rotational speed of blades 

Subscripts : 

f frame fixed with respect to undisturbed medium 

m moving frame fixed to aircraft 

r radiation direction 

ret evaluated at T = T* 

1,2,3 coordinate directions 


THEORY AND EQUATIONS 
Theoretical Formulation 

The problem of noise prediction can be reduced to the solution of the wave 
equation with a given distribution of sources on a moving boundary (the pro- 
peller blade surface) . There are basically two steps in the prediction process 
as follows: 

(1) The determination of the source distribution on the moving boundary 

(2) Solving the three-dimensional wave equation with specified source dis- 

tributions on the blade surface 

Previously formulated noise prediction theories (refs. 1 to 4) have been 
somewhat successful with conventional propeller designs, but all these methods 
have certain limitations which make them unsuitable for applications to pro- 
pellers with advanced geometry. Some of these limitations are the assumption 
of compactness of the acoustic sources and the far-field positioning of the 
observer. Three new requirements were considered in the development of the 
present prediction technique as follows: 

(1) Sources on the blade surface could not be approximated by sources in a 

plane because the blades have considerable twist and out-of-plane 
sweep 

(2) The propeller helical tip speed may be supersonic 

(3) Near-field noise calculations should be included 
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These requirements led to the selection of a purely numerical approach for 
noise calculations. The formulations used are two forms of the solution of the 
Ffowcs Williams-Hawkings (FW-H) wave equation (ref- 5) where the quadrupole 
term has been dropped- Thus, only the surface source terms have been retained - 
These solutions are derived in reference 6. 

^ The governing equation for the determination of the acoustic pressure 
p*(x,t), the Ffowcs Williams-Hawkings equation, is 


" 


V^P' 





( 1 ) 


where c and are the speed of sound and the density in the undisturbed 

medium, respectively. The variable v^ is_^the local normal velocity on the 
blade surface described by the equation f(x,t) = 0. The local force per unit 
area exerted by the blade surface on the fluid is denoted by and 6(f) is 

the Dirac delta function. 

The two forms of solution for equation (1) which are used to calculate 
propeller noise are 
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In equations (2) and (3) , dS is an element of surface area on the blade f = 0, 
and dT is an element of length of the curve of intersection of the surfaces 
f = 0 and g = T - t + (r/c) = 0- The symbol 0 represents the an^le between 
the normal to the surface f = 0 and the radiation vector r = x - y. In 

equation (2) , M = v * r/c, where v is the local velocity on the blade sur- 

face, r = r/r is the unit vector in the radiation direction, and Z-j- = ^j_rj_ 

is the force per unit area on the fluid in the radiation direction. The col- 
lapsing sphere method discussed in this paper refers to the application of 
equation (3). (See ref. 6.) 


The computer program discussed in this paper can handle cases in which the 
observer has the same_^forward speed as the propeller. This capability is 
achieved by letting Xj^ be the position of the observer in the moving frame; 
therefore, the observer position in the frame fixed to the undisturbed medium 
is 
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v(t*) dt* 



where v(t*) is the forward velocity of the propeller. Therefore, in the 
moving frame, the numerical value of p* (Xj^,t) is the same as p* (x^,t). In 
this instance, x^ = Xf{x^,t), and the value can be determined if the motion 
of the propeller is specified. 


Numerical Approach 

Figure 1 shows the flow chart for program PROPFAN. To commence the noise 
calculation, each blade is divided into panels. Since the blades of advanced 
propellers can have substantial amounts of twist and sweep, the noise sources 
may not be assumed to be in a plane. For this reason, a coordinate frame as 
shown in figure 2 is used to account for the three-dimensional character of the 
blades. This frame, called the rotating n-frame, is fixed to each blade. The 
center of this frame is on the propeller axis, and the coincides with 

the pitch change axis (PCA) of each blade. 

In order to describe the blade, the leading-edge curve of the blade is 
given as a function of distance along the pitch change axis r) 2 * Then the 
chord, airfoil section, twist, and thickness ratio of the blade are specified 
as functions of f )2 - The specification of these parameters is sufficient to 
describe the blade completely. 

To subdivide the blades into panels, each blade is cut in the radial direc- 
tion by planes perpendicular to the pitch change axis. In the chordwise direc- 
tion, a new nondimensional variable Q is introduced. This variable is the 
distance from the leading edge along the chord divided by the local chord CH, 

The upper and lower surfaces of the blade are now mathematically described as a 
function of the variables ri 2 and Q. From these values the surface position 
vector, unit normal vector, and surface area for each panel can be calculated. 
Figure 3 gives the definition of some of the parameters and symbols for the 
blade description. 

For the thickness noise calculations , both the upper and lower blade sur- 
faces are divided into panels. For loading noise calculations, the same blade 
division method is used if the upper and lower surface pressures are specified. 
However, if the pressure differential is specified, the mean surface (comprised 
of local chord lines) is divided into panels- This option between describing 
the input pressure function as either surface pressure or differential pressure 
is available to the program user because details of blade surface pressure dis- 
tribution are not generally available. 

After each blade has been divided into panels, the value of the helical 
Mach number of each panel is calculated to determine whether the panel is in 
subsonic or transonic/supersonic operation. Two different methods are used to 
determine the emission time. These two methods are explained in appendix A. 

If the element has subsonic motion, method 1 is used. If the panel moves at 
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transonic or supersonic tip speed, the possibility of multiple emission times 
exists and the more complicated method 2 must be used. 

The contributions of the ith panel to the acoustic pressure, denoted 
Pi * r from equations (2) and (3), may be written in finite-difference form 
as 


4np.'<J,t) = ifjl 
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ret 
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ret 
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Once the emission time (or times) of a panel has been calculated, the value 
of the radiation Mach number in the observer direction at the panel center is 
calculated- If this Mach number, M^, is less than a prescribed input value 
(usually taken as 0.98), equation (4) is used for that panel. For larger values 
of M^, equation (5) is used. The combined use of these two formulas improves 
the accuracy and efficiency of program execution. Although the application of 
equation (4) is simpler and faster than equation (5) , it cannot be used as the 
value of M^ approaches unity - 

Finally, the contributions of all the panels are summed to obtain the 
acoustic pressure signature p*(x,t). The numerical differentiations in equa- 
tions (4) and (5) are performed after the summation is completed . 


PR0GR?^M DESCRIPTION 

Program PROPFAN consists of a main program element and 19 additional sub- 
routines (including 3 input subroutines) . The program was written for batch 
execution but can be converted to interactive format if desired. The program 
input includes three subroutines and a FORTRAN namelist of input parameters. 

The program is written in FORTRAN IV and has been implemented at the 
Langley Research Center under the Network Operating System (NOS) on Control Data 
CYBER 170 series computer systems. The program as used at Langley has a maximum 
central memory requirement of approximately 43535 (octal) 60-bit words. A 
typical single-case run for a propeller operating at subsonic tip speed requires 
300 CPU seconds. For propellers with supersonic tip speeds, depending on the 
number of blade passage harmonics required, a typical case may take several 
thousand CPU seconds. 

There are two versions of the program currently in use. The version called 
PROPFAN developed for use by industry has only printed output capability. In 
view of the great diversity in computing and graphics hardware, it was felt that 
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this basic version of the program will be more valuable for easy implementation 
on users* computing systems- Users may add their own routines to program 
PROPFAN for plotting, data storage, and editing options according to their own 
needs . 

The version presently used at Langley is called NPROPFN and has elaborate 
graphics capability. An example of a control card deck used to execute this 
version at Langley is given in appendix B. Examples of the printed output of 
PROPFAN are presented in appendix C. The graphic output of NPROPFN for the 
examples given in appendix C are presented in appendix D. The two versions 
of the program are mathematically identical in the main part of the code where 
the acoustic calculations are performed. 

The following sections describe each program element and its function for 
PROPFAN. Note that for identification purposes only, the name of each subrou- 
tine in NPROPFN is kept the same as that in PROPFAN except for the letter N 
preceding the first letter. For example, the subroutine KRUD in PROPFAN is 
called NKRUD in NPROPFN. 


Main Program - PROPFAN 

PROPFAN , the main program, is the driving element which calls each main 
subroutine in the proper order. Calls INTLZ, GMTRY, DBYDT, and DOPRINT. 


Main Subroutines 

INTLZ initializes the input parameter default values and reads in namelist 
input changes. Some constants are calculated and storage arrays are zeroed. 
Called by PROPFAN. 

GMTRY divides each blade into panels using spanwise and chordwise divisions 
according to the input values of NPCA, NPCAF, NLE, NCHF, and THKMAXL. The geom- 
etry and coordinates of each panel are calculated after calling the input sub- 
routines FUNE2 and FUNE2Q- The blade pressure functions are obtained from input 
subroutine FUNPRES - Each panel's contribution to the acoustic pressure at each 
observer time is found by calling PROP. The torque, thrust, and power are also 
calculated in GMTRY. Called by PROPFAN; calls FUNE2, FUNE2Q, FUNPRES, and PROP. 

PROP steps each panel through observer time and adds each contribution to 
the total acoustic pressure. Three pressure arrays are generated: thickness, 

loading, and drag (skin friction). The near- and far-field terms are calculated 
for loading and drag noise. Subroutines QINTLZ and QULOCK are first called for 
an initial estimate of the time step. Called by GMTRY; calls QINTLZ, QULOCK, 
QURUD, KRUD, and KLOCK. 

DBYDT evaluates the integrals in the acoustic equations and performs dif- 
ferentiation with respect to the observer time if required (refer to eqs . (4) 

and (5) ) . This routine produces the final output consisting of the component 
pressure signatures and spectra. It also contains code which rotates the 
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acoustic pressure signature (if periodic) to center the pressure peak. This is 
useful if graphics output is used. The acoustic pressure spectra are calculated 
in this subroutine. Called by PROPFAN. 

DQPRINT generates the printed output- The basic output consists of a Data 
Sheet listing the input parameters and output values of torque, thrust, power, 
and sound pressure level. A more detailed output is optionally available (con- 
trolled by logical variable PRINT) which includes the pressure signatures and 
spectra of the thickness, loading, and drag noise components. Called by 
PROPFAN. 

ADDIN adds separately the positive and negative contributions of the inte- 
grals in equations (4) and (5) to improve the accuracy of the integrations. 
Called by CORRECT, KRUD, and QURUD. 


Subsonic Logic 

KRUD calculates the pressure integrands and sums them for blade panels 
which move at subsonic speed. Called by PROP; calls KLOCK and ADDIN. 

KLOCK calculates the retarded time for any observer time for blade panels 
which move at subsonic speed. Called by KRUD. 


Supersonic Logic 

QINTLZ initializes the parameters for the first call to QULOCK. Called by 

PROP. 


QURUD calculates the pressure integrands and sums them for blade panels 
which move at transonic/supersonic speed. Called by PROP; calls QULOCK and 
ADDIN. 

QULOCK calculates the retarded time or times for a fixed observer time for 
blade panels moving at transonic/supersonic speed. Called by QURUD. 

CORRECT calculates the pressure integrands and sums them for blade panels 
which move at transonic/supersonic speed- This part of the code corrects for 
the l/|l-M^| singularity (eq. (4)) as Mj. approaches unity by replacing it 
with an equivalent expression using the collapsing sphere technique. Called by 
QURUD; calls QINTLZ, QULOCK, GETSET, CROSSED, and ADDIN. 

CROSSED determines if the collapsing sphere centered at the observer 
crosses the blade panel for supersonic regions' of the blade; if so, the width 
of the crossing is calculated. Called by CORRECT. 

GETSET calculates the coordinates of the corners of each blade panel for 
use in subsequent calls to CROSSED. Called by CORRECT; calls FUNE2. 
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Interpolation Routines 


PCS perforins curve interpolation from a data array using the piecewise 
cubic spline method. Called by FUNE2, FUNE2Q, or FUNPRES if input functions are 
in tabular data form; calls ORTHO. 

ORTHO calculates the coefficients of an Nth order polynomial which provides 
the best least-squares fit for a set of points. Called by PCS. 


OPERATING INSTRUCTIONS 

The program PROPFAN requires certain input functions and a namelist input. 
In the following sections these requirements are described. 


Input Functions 

Three input subroutines, FUNE2, FUNE2Q, and FUNPRES, are required by pro- 
gram PROPFAN. These subroutines are used to describe the physical and aerody- 
namic characteristics of the propeller blade . 

Subroutine FUNE2 defines the blade angle of attack, leading-edge curve, 
chord width, and thickness ratio as functions of the location along the blade 
span only. Routine FUNE2Q defines the camber and thickness functions of the 
blade as each varies with location along the span and/or chord. The third, 
FUNPRES, describes the blade pressure distribution functions which may also 
depend on span and/or chord coordinates. The pressure function may be described 
as surface or differential pressure as is convenient to the user. 

These three routines allow great flexibility in the definition of blade 
geometry and loading. The input variables may be described by using analytic 
functions or they may be input in tabular data form. If the input is in array 
form, a mathematical interpolation routine (such as routine PCS) must be used to 
calculate the required set of points. 

Subroutine FUNE2 : The variables in this subroutine are all functions of 

radial distance E2 along the pitch change axis. Note that R is the propeller 
radius. Calling sequence: (E2, R, AA, AAP, CH, CHP, LED, LEDP, THKRAT, 

THKRATP, CAMRAT, CAMRATP) . Refer to figures 2 and 3. 

Variable Definition 


AA Real. Angle of attack (degrees). Angle between mean 

chord line and Uj^-axis . 

CH Real. Local chord (meters). Distance along mean chord 

line from leading edge to trailing edge. 
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LED 


Real. Leading-edge displacement (meters). Distance from 
leading edge to pitch change axis, measured along line 
parallel to mean chord line. Positive direction of 
measurement is from leading edge to trailing edge. LED 
is typically negative. For conventional straight blades, 
use LED = -0.25*CH. 


THKRAT Real. Airfoil thickness ratio. Ratio of maximum airfoil 

thickness at each radial location to local chord . 

CAMRAT Real. Airfoil camber ratio. Ratio of maximum camber at 

each radial location to local chord. 


AAP, CHP, LEDP, 
THKRATP, CAMRATP 


Real. Derivatives of the above variables with respect to 

8 1 3 

E2. Note that where ER = E2/R (nondimen- 

sional radial distance) . 


Subroutine FUNE2Q : The variables in this subroutine are generally func- 

tions of nondimensional chordwise variable Q. The variable E2 (radial dis- 
tance, meters) is carried into this subroutine to allow for any change of air- 
foil section at different radial stations. Variables in this subroutine can 
also be functions of E2 if required. Calling sequence: (E2, R, Q, CMBR, CMBRP, 

THK, THKP) . Refer to figure 3. 


Variable 

CMBR 


THK 


CMBRP, THKP 


Definition 


Real. Nondimensional camber. Distance from chord surface 
to camber surface, divided by camber ratio times chord. 
Generally function of Q, where Q = X/CH. For symmetric 
airfoil, CMBR equals zero. 

Real. Nondimensional thickness. Distance from camber 
surface to blade surface, divided by thickness ratio 
times chord. Generally a function of Q. 

Real. Derivatives of above variables with respect to Q. 


Subroutine FUNPRES : This routine describes the pressure and viscous shear 

stress distribution on the blade surfaces. These distributions can be a func- 
tion of the radial distance E2 and/or the nondimensional chord location Q. 


The input may be expressed as pressure on the upper and lower surfaces 
(absolute pressure minus atmospheric pressure) or as a differential pressure. 

If the surface pressure approach is used, the user can include a description of 
the viscous shear stress (skin friction) for use in the noise calculations. If 
the input is in differential pressure form, the program does not include drag 
terms in the calculations and these variables must be set equal to zero. 

Calling sequence: (E2, R, Q, SURFP, DP, SPU, SPL, SIGMAU, SIGMAL) . 
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Variable 


Definition 


SURFP Logical. If true, the variables being returned to the pro- 

gram are upper and lower surface pressures (SPU and SPL) 
and skin friction stress (SIGMAU and SIGMAL) may be 
defined or set equal to zero. If false, the variable 
returned is pressure differential (DP) and skin friction 
stress values must be set equal to zero. 

DP Real. Differential pressure (pascals) acting on camber 

surface. This variable may be defined as DP = SPL - SPU. 


SPU, SPL Real. Absolute pressure minus ambient pressure acting on 

upper and lower blade surfaces, respectively (pascals). 

SIGMAU, SIGMAL Real. Skin friction (viscous shear) stress acting on upper 

and lower blade surface (pascals) , Set equal to zero if 
SURFP is false. 


Namelist Parameters 

A namelist called $INPUT is used to read in the basic parameters such as 
blade dimensions, rotational speed, and forward propeller speed. The namelist 
also includes several parameters which control the speed and accuracy of program 
execution. For example, the number of spanwise and chordwise divisions, number 
of harmonics, the allowable time step error, and similar variables are defined 
through the namelist. In addition, there are four logical flags which control 
the specific amount and type of printed output. 

For each namelist variable, there is a default value defined within the 
program which applies to the case of a three-blade propeller in subsonic opera- 
tion. If any of these variables are not redefined by the namelist input during 
a particular single-case run, the default values are automatically used. 
Multiple-case runs which vary the namelist parameters while holding the input 
subroutines constant can be executed within the same job. This may be done by 
submitting sequential namelists and identifiers for each case. However, once a 
variable is defined within a job, that value will be retained for all the fol- 
lowing cases within that job until it is subsequently redefined. 

A complete list of the namelist parameters, their definitions, and default 
values is given in the following list. An example of the control cards and 
input format for a multiple-case run of program NPROPFN (including graphics 
package) may be found in appendix B. 
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Variable 

C 

CHECK 


CORCALL 

EPSILON 


FRACDT 

MOTION 

MTRANS 


NBLADES 


Definition 


Real. Speed of sound. 

Logical. When true, the following intermediate 
values from subroutine GMTRY are included as 
output: iteration indices for spanwise and 

chordwise stepping; geometry of each blade 
element; coordinates of each element 

on the mean, upper, and lower surfaces; unit 
normal vector components on each surface. 

The output is particularly useful for 
debugging . 

Logical. When true, the number of calls to 
CORRECT for each observer time is included 
in printed output. 

Real- Maximum allowable error in iterative 
calculation of retarded time, defined as a 
percent of observer time step used. The 
iteration process for finding retarded time 
will stop when two successive approximations 
of retarded time differ by less than 
0-01 X EPSILON X At where At is observer 
time step in the program. Reducing 
EPSILON will increase execution time. 

Real. Fraction of time step to use in 
observer time differentiation- For sub- 
sonic cases, use 0.5 = FRACDT = 1.0; for 
supersonic blades, use 1.0 = FRACDT = 3.0. 

Logical. When true, observer moves forward 
with propeller. When false, observer is 
stationary with respect to undisturbed 
medium - 

Real. If helical Mach number of blade panel 
becomes greater than or equal to MTRANS, 
control is switched from subsonic to 
transonic/supersonic logic . For efficient 
program execution, use a range of 0.95 
to 0.99. 

Integer. Number of propeller blades. 


Default 
345.0 m/sec 
False 


False 

0.5 percent 


0.50 

True 

0-95 
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2 


NCHF 


NLE 


NPCA 


NPCAF 


NPTS 


NSPEC 


NTAU 


NTE 


PLOTS 


Integer. Factor to multiply chordwise divi- 
sions (NLE and NTE) for transonic/supersonic 
regions of blades in order to create a larger 
number of smaller panels. 2 = NCHF = 4 is 
recommended. 

Integer. Number of chordwise divisions from 
leading edge to THKMAXL. These divisions 
are unequally spaced by code to force a 
greater number of divisions close to lead- 
ing edge. For sharp leading edge, set 
NLE = 0. For blunt leading edge, 

5 = NLE = 10 is recommended. 

Integer. Number of spanwise divisions for 
subsonic regions of blade. Optimum range 
is 15 = NPCA = 25. 

Integer. Factor to multiply number of span- 
wise divisions (NPCA) for transonic/ 
supersonic regions of blades in order to 
create larger number of smaller panels. 

2 = NPCAF ^ 4 is recommended. 

Integer- Number of points in time per blade 
passage period. For subsonic blades, 

100 = NPTS = 200 is adequate. For super- 
sonic tip speeds, use rule NBLADES*NPTS<650 . 
The value of NPTS must always be even and 
less than 400. 

Integer. Number of harmonics used to calcu- 
late acoustic spectrum. NSPEC must be less 
than one-half of NPTS. 

Integer. Number of source time subdivisions 
to use in collapsing sphere routine for 
panels moving at transonic/supersonic 
speed. Increasing NTAU increases accuracy 
of acoustic calculations but also increases 
program execution time. 

Integer. Number of chordwise divisions from 
THKMAXL to trailing edge. These divisions 
are equally spaced by the code. A good 
range for this variable is 10 = NTE = 20, 

Logical. Used only in NPROPFN version. When 
true, plots of output data are generated. 
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150 


30 


10 


10 


False 
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PRINT Logical. When false, only Data Sheet is True 

printed. When true, component pressure 
signatures and spectra are also printed. 


R 

Real . 

Blade radius . 

1.3 m 


REV 

Real . 

Propeller rotational speed - 

2145.0 

rpm 

RHO 

Real. 

Density of air- 

1.2029 

kg/m^ 

RINNER 

Real - 

inner blade radius. 

0.3 m 


THKMAXL 

Real . 

Location of transition point between 

0.05 



NLE and NTE type chord divisions, defined 
as fraction of chord from leading edge. 
For a sharp leading edge, set THKMAXL = 0 
and NLE = 0. For other cases, use 
0.05 = THKMAXL = 0.10. 


TRANS 


TWICE 

V3 


XO 


Real. If radiation Mach number (M^^) from an 0-98 

element to observer becomes greater than 
TRANS, noise calculations are done using 
collapsing sphere technique (eq. (5)). 

Logical. When true, output is printed twice. False 


Real- Propeller forward speed, perpendicular 40.0 m/sec 

to plane of rotation (xiX 2 -plane) . 

Real, three dimensional. For moving observer, (7.28,0.,0.) m 

this is initial position with respect to 
frame fixed to propeller center. Propeller 
plane is in xiX 2 “plane and axis of propeller 
is along X 3 -axis. For stationary observer, 
observer must be specified in frame which 
has origin at center of propeller, with same 
orientation, at observer time t = 0. This 
frame will remain stationary for later 
observer times . 


Output Identification 

A card with a 10-character identifier (Format AlO) in the first 10 columns 
is used following the namelist $INPUT in order to label each case run. Failure 
to use this identifier will produce an error message. The output Data Sheet, 
pressure signature, and spectrum tables will all be labeled with this 
identifier. 
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Output Data 


The output data initially consists of a Data Sheet listing each of the 
namelist input parameters and the calculated output values of torque, thrust, 
power, and component sound pressure levels. The Data Sheet also lists the 
value of the observer time step (DT) for the calculation of acoustic pressure 
and the number of calls (INDEXSL) to the supersonic code (subroutine CORRECT) . 

Two tables listing the pressure signature (in pascals) and the pressure spectrum 
(in decibels) are optionally available through the use of the logical input 
variable PRINT. 

The pressure signature output is the time history of the thickness, loading, 
drag, and overall noise levels calculated for the specified number of observer 
times (input variable NPTS) . When the acoustic pressure is periodic, the pro- 
gram rotates the signature to center the pressure peak within the data array. 

This rotation is useful in producing well laid-out plots of the output. At the 
bottom of the signature table, one point number is labeled with the letter C. 

This indicates that this point number (with the corresponding pressure values) 
was actually the initial zero observer time used to start the noise calcula- 
tions before rotation of the data array. Thus, relative phase information can 
be obtained by comparing the starting points of different case runs. If the 
observer is stationary and the propeller is in forward motion, the pressure 
signatures are not periodic. 

The pressure spectrum output lists a table of the component and overall 
noise levels calculated for the specified number of harmonic frequencies (input 
variable NSPEC) . 


CONCLUDING REMARKS 

The computer program discussed in this report is based on two of many 
theoretical formulations available for prediction of the noise of rotating 
blades. These formulations were selected for their simplicity and ease of pro- 
graming. The blade description in the program is such that complicated designs 
for advanced propellers can be specified easily. All geometric parameters of 
the blade such as the twist and spanwise variation of blade thickness can be 
prescribed precisely and utilized in the acoustic calculations. Comparisons of 
the output of this program with measured acoustic data have been very encourag- 
ing. It is hoped that this program will be beneficial to the users as a tool 
in propeller design and noise study. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
June 12, 1981 
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Figure 1.- Flow chart of propeller noise prediction 
program PROPFAN. 
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E2 

radial location on ri2“9xis 


AA 

geometric angle of attack (functi 

on of TI 2 ) 

PCA 

pitch change axis 


CH 

local chord (function of 02 ) 


Q 

distance from LE divided by local 

chord (CH) 


Figure 2.- Curvilinear coordinate system (Qf^2^ used to describe 
blade geometry of advanced propellers. Blade mean surface is 
not in plane. 
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CALCULATION OF EMISSION TIME 

If the source and the observer times are T and t, respectively, and the 
distance between the observer and the source is r, the emission time T = T* 
is calculated from the relation 


c(T* - t) + r = 0 


(Al) 


The emission time is the time when the sources on the panel emit sound which 
arrives at the observer at time t. The source position y is a function of 
T, so that r itself is a function of T. Because of the trigonometric terms 
in r (eq. (A4) ) , equation (Al) cannot be solved for T* in closed form. A 
numerical method must be used. Two schemes are used depending on the speed of 
the source . 


Method 1: Subsonic Helical Tip Speeds 

When the observer time and position are fixed, one can show for a given 
source in motion that 

If = 1 - Mr (A2) 

where g = x - t + (ir/c) . Equation (Al) indicates that, viewed as a function of 
a single variable T, the emission times of a source in motion are the zeros of 
function g. For sources in subsonic motion, one has Mj. < 1 and 8g/9x > 0. 
This means that the function g is a strictly increasing function of T and 
thus can have only one zero. 

Newton's method was used to develop a fast iterative scheme to find the 
emission time of each panel on the blade. To speed up the convergence of 
iteration, the known emission time of a nearby panel is used as the initial 
guess for the emission time of the next panel. This method has resulted in 
shortened computation time for general aviation propellers operating at sub- 
sonic tip speeds . 


Method 2: Transonic/Supersonic Helical Tip Speeds 

For panels moving at supersonic speeds, the range of search for each root 
of equation (Al) must be narrowed to find the appropriate nontrivial solutions. 

Assume that the source lies in the xiX 2 “plane at source time T = 0 . The 
source position at this time is described by (H/^) in polar coordinates as 
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shown in figure Al. The source rotates around the X 3 ~axis (propeller axis) at 
an angular velocity of oj and moves forward in the positive direction of the 
axis at velocity V 3 . The x-frame is fixed to the undisturbed medium and is not 
in motion. The observer position and time are x and t, respectively. 

Viewed by the observer in this frame, the source path is a helix. From equa- 
tion (Al) , the following equation may be obtained: 

c^(T* - t)2 - = 0 (A3) 

with only the solutions T* = t being physically acceptable. The square of the 
distance between the source and the observer can be written as follows: 


r^ = X 3 ^ + d^ + 2 dr) cos ((j) - 3) + V 3 ^X^ - 2 X 3 V 3 T (A4) 

where (}) = + OJT, and d is the distance of the observer from the X 3 ~axis. 

(Refer to fig. Al . ) Assiaming V 3 /C = M < 1 and Tq = X 3 /C, equation (A3) can 
now be written as 

K(T) = c^(l - M^)T^ - 2c^(t - MTq)T + C^(t^ - Tq^) 

= d^ + “ 2 dn cos (4) - 3) (a5) 


The graphs of 

K(T) = C^(l - M^)T^ - 2c^(t - MTq) + C^ (t^ - Tq^) 
K(t) = + r|^ - 2 dr| cos ((|) - 6 ) 


are a parabola and a sinusoidal curve in the TK-plane, respectively. Thus, to 
find the roots of equation (A5) , the points of intersection of a parabola and a 
sinusoidal curve must be determined. These roots must lie between the points 
of intersection of the lines K = (d + p) parallel to the T-axis in the 
TK-plane. These lines are obtained by setting cos (({) - 3) in equation (A5) 
equal to -1 and 1 , respectively. 

Not all the roots of equation (A5) are emission times because the require- 
ment T* = t is not satisfied for all the roots. Only the roots obtained by 
the intersection of the parabola and the sinusoidal curve which are to the left 
of the axis of the parabola are the emission times. One can graphically demon- 
strate that only an odd number of emission times are obtained if multiple roots 
are counted by their multiplicity (e.g., two equal roots are counted as two 
roots) . The method of modified regula falsi discussed in reference 7 is used to 
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find the roots of equation (A5) . In the implementation of this method ^ the 
intervals in which each root of the equation lie are found first. An iterative 
method is then used in each interval to find the roots. This method is more 
time consuming to implement on the computer than that of method 1. However, of 
the several methods for finding emission time tried in the course of the 
development of the computer program, the scheme based on modified regula falsi 
was the most satisfactory, both in precision and speed and in guaranteeing that 
all the roots would be found. 
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EXAMPLE CONTROL CARD DECK 

This appendix presents an example of a control card deck used at Langley 
Research Center to execute a multiple-case run of program NPROPFN including the 
graphics package. 


* JOBNAME, ACCOUNT AND CHARGE CARDS 

* 


JOBNAME, T50000. 

ACCOUNT, 123456N. 

CHARGE, 987654, LRC. 

* 

* COPY INPUT SUBROUTINES TO APPROPRIATE FILES, 

* COPY NAMELIST INPUT TO TAPE! . 

* 

COPYBR, INPUT, FUNE2. 

C0PYBR,INPUT,FUNE2Q. 

COPYBR, INPUT, FUNPRES. 

REWIND , FUNE2 , FUNE2Q , FUNPRES . 

GET,ND0PR0P/ST=CPF,UN=564526N. 

CALL,ND0PR0P. 

* 

* END OF RECORD, READ IN FUNE2. 

* 

7-8-9 

SUBROUTINE FUNE2 

$ ( E2 ,R ,AA ,AAP ,CH ,CHP ,LED ,LEDP ,THKRAT ,THKRATP ,CAMRAT ,CAMRATP ) 

I I 


RETURN 

END 

* 

* END OF RECORD, READ IN FUNE2Q. 

* 

7-8-9 

SUBROUTINE FUNE2Q 
$ (E2,R,Q,CMBR,CMBRP,THK,THKP) 

I I 


RETURN 

END 

* 

* END OF RECORD, READ IN FUNPRES. 
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SUBROUTINE FUNPRES 

$ (E2,R,Q,SURFP,DP,SPU,SPL,SIGMAU,SIGMAL) 

I I 


RETURN 

END 

* 

* END OF RECORD, READ IN NAMELIST INPUT. 

* 

7-8-9 
$ INPUT 
C = 344.7, 

NLE = 8, 

NPTS = 200, 

NTE = 12, 

RHO = 1.206, 

V3 = 63.79, 

XO = 7.28, 0., 0. $ 

EXAMPLE 1 
$ INPUT 

XO = 7.28, 0., -3.05 $ 

EXAMPLE 2 
* 

* END OF INFORMATION. 

* 

6-7-8-9 
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NDOPROP PROCEDURE LISTING 


* START OF NDOPROP 

* (INPUT ON TAPEl, OUTPUT ON TAPE2) 

* 

* GET MAIN PROGRAM AND READ INPUT SUBROUTINES. 

* 

GET , NPROP FN/ST=CP F , UN=564526N . 

COPYBR,NPROPFN,LIST. 

C0PYBR,FUNE2,LIST. 

C0PYBR,FUNE2Q,LIST. 

COPYBR,FUNPRES,LIST. 

PACK, LIST. 

* 

* COMPILE MAIN PROGRAM AND INPUT SUBROUTINES 

* 

FTN,I=LIST,L=MYLIST,B=BINARY,A,P. 

GOTO, 200. 

EXIT. 

REWIND, MYLIST. 

COPY, MYLIST, OUTPUT. 

EXIT. 

* 

* EXTRACT INPUT FUNCTIONS FROM FUNE2,FUNE2Q,FUNPRES 

* 

200 ,REWIND,FUNE2 ,FUNE2Q ,FUNPRES. 
C0PYBR,FUNPRES,TAPE4. 

C0PYBR,FUNE2Q,TAPE4. 

C0PYBR,FUNE2,TAPE4. 

PACK TAPE4. 

GET , FUNEX/ST=CPF ,UN=564526N . 

EDIT,TAPE4,N,FUNEX,0. 

RETURN, NPR0PFN,FUNE2,FUNE2Q,FUNPRES, LIST. 

* 

* ATTACH LIBRARIES, LOAD AND EXECUTE 

* 

GET,NPR0LIB/ST=CPF,UN=564526N. 

GET,MYLIB4/ST=CPF,UN=564526N. 

ATTACH , F I NML I B , LRCGOS F/UN=L I BRARY , NA . 
LDSET,LIB=NPR0LIB/MYLIB4/FTNMLIB/LRCG0SF. 

LDSET ,MAP=SB/PROPMAP , PRESETA=DEBUG . 

LOAD, BINARY. 

NOGO.FILE. 

RETURN , B I NARY , NPROL I B , MYL I B4 , FTNML I B, LRCGOS F. 

FILE. 

GOTO, 300. 

EXIT. 

REWIND, PROPMAP. 

COPY, PROPMAP, OUTPUT. 

EXIT. 
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* RETURN FILES 

* 

300, RETURN, TAPES, TAPE?, FILE. 

* 

* PLOTTING INSTRUCTIONS. 

* 

PLOT.CALPOST.FB (FRAMCNT=1 20) 

CONT. // SHEET PAPER - 8-1/2 x H INCHES 
CONT. BLACK INK. LEROY .3 MM 

CONT. 

CONT. X=8.5 AND Y=ll // 

GOTO, 400. 

EXIT. 

* 

* END OF NDOPROP 
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EXAMPLES OF PRINTED OUTPUT 

In this appendix, four examples of the output of PROPFAN are presented. 
Examples 1 and 2 are for a conventional propeller in subsonic operation. 

Examples 3 and 4 are for an advanced design propeller in supersonic operation. 

A list of the input functions and parameters (Data Sheet) is always printed in 
the output listing. The noise component pressure signature and pressure spectra 
output which are optionally available are also included here. These examples 
were chosen to demonstrate the versatility of PROPFAN with regard to the input 
and output data. The graphic output (generated by NPROPFN) for these examples 
is presented in appendix D. Several other examples for advanced design pro- 
pellers and comparisons with measured noise data are given in reference 8. 


Examples 1 and 2 

Examples 1 and 2 illustrate the use of analytic functions to describe a 
conventional propeller. In subroutine FUNE2, the variables AA, CH, LED, and 
THKRAT are all given as functions of E2, while the camber ratio (CAMRAT) is 
defined as a constant. In subroutine FUNE2Q, the camber (CMBR) is defined as a 
function of the nondimensional distance along the chord line (Q) . In addition, 
two different thickness functions (THK) are specified which vary with Q and 
with radial location E2. 

In example 1, because logical variable SURFP is set false, subroutine 
FUNPRES supplies pressure as the differential pressure (DP) , and the skin 
friction values (SIGMAU and SIGNAL) are necessarily set equal to zero (refer to 
section "Input Functions," Subroutine FUNPRES). In example 2, logical variable 
SURFP is set true, and the subroutine returns the surface pressures (SPU and 
SPL) . In this case, SIGMAU and SIGNAL are defined as functions of radial dis- 
tance E2. The observer locations (XO) are also different in these two examples. 


Examples 3 and 4 

Examples 3 and 4 illustrate the use of arrays of input data to describe 
the geometry and loading of an advanced design propeller operating with super- 
sonic helical tip Mach number. In subroutine FUNE2, the data arrays and asso- 
ciated variables are labeled as follows: 

AAl angle of attack (degrees) , data array corresponding to ERl 

ERl radial location divided by propeller radius (nondimen- 

sionalized radial position) 

OFFSET offset (degrees) for angle of attack to obtain desired 

pitch angle at operating condition 
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CH2 

THKRT2 


chord divided by propeller diameter, data array corre- 
sponding to ER2 

thickness radio, data array corresponding to ER2 


LED 2 


leading-edge displacement divided by chord, data array 
corresponding to ER2 


ER2 


radial location divided by propeller radius (nondimen- 
sionalized radial position) 


WORKl, W0RK2, storage arrays for interpolation subroutine PCS 
WORKS, W0RK4 


II, 12, IS, I4 initialization indices for interpolation subroutine PCS 


The specific arrays needed for program execution are interpolated from these 
data arrays using subroutine PCS. Note that the array of LED2 values has the 
opposite sign as LED, because this is how the blade is specified for manu- 
facture. The derivatives AAP, CHP, and THKRATP are taken with respect to radial 
distance (ERl for AAP and ER2 for the remaining derivatives) . The variable 
CAMRAT is defined analytically as a function of radial distance E2 and the 
derivative CAMRATP is taken with respect to E2, 

In subroutine FUNE2Q, the airfoil thickness (THK) and camber (CMBR) func- 
tions are described analytically. Variable CMBR is defined as a function of Q 
and of E2 because the inner section of the blade (ER = 0.45) uses a series 65 
airfoil and the outer section (ER > 0.45) uses a series 16 airfoil section. 


In subroutine FUNPRES, the differential pressure (DP) across the mean sur- 
face of the blades is used for both examples. Thus, variable SURFP must be 
defined as false, and the skin friction stresses (SIGMAU and SIGMAL) must be 
set equal to zero. The pressure distribution (DP) is derived from a data array 
of lift coefficients. The additional variables are defined as follows: 


CP2 

ER2 

INTRVL 

FACTOR 

OMEGA 

RHO 

VINES 


lift coefficient, data array corresponding to ER2 

radial location divided by propeller radius (nondimen- 
sionalized radial position) 

initialization value for interpretation subroutine PCS 

scaling factor to adjust lift coefficient to obtain 
desired power absorbed by the blades 

propeller speed (radians/sec) 

density of air (kg/m ) 

relative wind velocity at blade tip, squared (m^/sec^) 


30 



APPENDIX C 


The lift coefficients at each radial location are interpolated from the CP2 data 
array using subroutine PCS. Then the differential pressure (DP) values are cal- 
culated by using the lift coefficients and nondimensional chord location Q. 
Thus, in this case the pressure varies with both radial and chordwise location. 

The only difference between examples 3 and 4 is the position of the 
observer (XO) . In example 3, the observer is located almost in the plane of 
the propeller and in motion with the propeller. Note that the supersonic sec- 
tion of the program is used many times during execution. This can be seen from 
the value of INDEXSL printed on the Data Sheet which indicates the number of 
calls to the supersonic code. In example 4, the observer is located behind the 
propeller and the blades do not move at supersonic speeds toward the observer 
at any time. In this case, the program executes very quickly since only the 
subsonic code is used. This efficiency is achieved mainly by using only coarse 
panel sizes (input variables NCHF and NPCAF are set equal to one) . 
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C 


c 

c 

c 


c 


c 


c 


c 


c 


c 

c - 
c 


c 


c 

100 


c 

200 


SUBROUTINE FUNE2 

$ ( E2 , R , AA , AAP , CH , CHP , LED , LEDP , THKRAT , THKRATP , CAMRAT , CAMRATP ) 

REAL LED, LEDP 
ER = E2/R 

TWIN OTTER INPUT FUNCTIONS, EXAMPLE 1. 

AAEXP = EXP(-3.491*ER) 

AA = 12.622 + 78.037*AAEXP 
AAP = -209.529*AAEXP 

ROOT = SQRT( 1.-ER ) 

POLY = 1.474 + ER»( 2.544 + ER ) 

CH = 0.069»ROOT»POLY 

CHP = O.069/R*(-0.5/ROOT»POLY + ROOT»(2.»ER+2.544) ) 

LED = -0.25»CH 
LEDP = -0.25*CHP 

THKEXP = EXP(-11.200»ER) 

THKRAT = 0.069 + 3.2244»THKEXP 
THKRATP = -27.778»THKEXP 

CAMRAT =0.03 
CAMRATP = 0. 

RETURN 

END 


SUBROUTINE FUNE2Q 

$ ( E2, R, Q, CMBR, CMBRP, THK, THKP ) 

TWIN OTTER INPUT FUNCTIONS, EXAMPLE 1. 

CMBR = 4.»(Q-Q»»2) 

CMBRP = 4.-8.»Q 

IF (Q.GT.0.5) GO TO 100 

THK = 0.989665*SQRT(Q) + Q»(-0. 23925 + Q»(-0.041 - Q»0.5594)) 
THKP = 0.4948325/SQRT(Q) - .23925 + Q»(-0.082 - Q»1.6782) 

GO TO 200 

QM = 1. - Q 
QMP = -1. 

THK = 0.01 + QM»(2.325 + QM»(-3.42 + QM»1.46)) 

THKP = (2.325 + QM»(-6.84 + QM»4.38) )»QMP 

RETURN 

END 


EXAMPLE 1 


32 



APPENDIX C 


SUBROUTINE FUNPRES EXAMPLE 1 

$ (E2, R, Q, SURFP, DP, SPU, SPL, SIGMAU, SIGMAL) 

C 

LOGICAL SURFP 
C 

C - TWIN OTTER INPUT FUNCTIONS, EXAMPLE 1. 

C - DIFFERENTIAL PRESSURE. 

C 

SURFP = .FALSE. 

A = . 719618*10. »»6. 

C 

ER = E2/R 

IF (Q-0.05) 10,20,20 
10 SPQU = EXP(800.»Q»(Q-.1)) - 0.5 

SPQL r 0.9*EXP(500.»Q»(Q-.1)) - 0.4 
GO TO 100 
C 

20 IF ((Q-.05)»(Q-.5)) 30,30,40 

30 SPQL = -0.1421 

SPQU = -0.3647 
GO TO 100 
C 

40 SPQL = .2842*(Q-1.0) 

SPQU = .7293*(Q-1.0) 

C 

100 SPU = SPQU»ER»»2.833*SQRT(1.-ER)»A 
SPL = SPQL»ER»»2.833*SQRT(1 .-ER')»A 
C 

DP = SPL - SPU 
C 

SIGMAU = 0. 

SIGMAL = 0. 

C 

RETURN 

END 
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DATA SHEET 


c = 

344.7 

M/SEC 




RHO = 

1.2060 KG/M»»3 




XO = 

( 7.2fi 

1, 0.00, 0.00) M OBSERVER IS IN MOTION 

V3 = 

63.8 

M/SEC 




= 

142.7 

MPH 




R = 

1.300 

M 

HELICAL M-TIP = 

.867 

RINNER = 

.300 

M 




REV = 

2145.0 

RPM 

PER = 

9.3240E-03 

SEC 

NBLADES = 

3 



BPF = 107.: 

3 HZ 

NPCA = 

15 



NPCAF 

= 3 

NLE = 

8 



NCHF 

= 2 

NTE = 

12 





THKMAXL = 

.05 



NTAU 

= 10 

NPTS = 

200 


DT = 

4.6620E-05 

SEC 

NS PEC = 

30 



FRACDT = 

.500 





MTRANS = 

.950 

EPSILON = 

.50 

% 


TRANS = 

.980 





INDEXSL = 

0 

TORQUE = 

581.43 

N-M/BLADE 




THRUST = 

1947.37 

N/BLADE 




POWER = 

130.60 

KW/BLADE 




= 

175.07 

HP/BLADE 




OASPL = 

113.52 

DB RE 20.E-6 PA 




= 

108.74 

DB RE 20.E-6 PA 

(THICKNESS) 



= 

111.69 

DB RE 20.E-6 PA 

(LOADING) 



= 

-158.57 

DB RE 20.E-6 PA 

(DRAG) 




EXAMPLE 1 
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PRESSURE SIGNATURES OF 


POINT 

OBSERVER 

THICKNESS 

— 

LOADING NOISE 

NUMBER 

TIME 

NOISE 

FAR FIELD 

NEAR FIELD 


(MSEC) 

CPA) 

(PA) 

(PA) 

1 

0.00 

.9814 

1.7755 

-2.1219 

2 

.05 

.9812 

1.8787 

-2.1214 

3 

.09 

.9813 

1.9819 

-2.1206 

4 

.14 

.9818 

2.0850 

-2.1196 

5 

.19 

.9826 

2.1882 

-2.1183 

6 

.23 

.9837 

2.2914 

-2.1168 

7 

.28 

.9853 

2.3947 

-2.1151 

8 

.33 

.9871 

2.4982 

-2.1132 

9 

.37 

.9894 

2.6019 

-2.1110 

10 

.42 

.9920 

2.7059 

-2.1087 

11 

.47 

.9950 

2.8101 

-2.1061 

12 

.51 

.9983 

2.9147 

-2.1032 

13 

.56 

1.0021 

3.0196 

-2.1002 

14 

.61 

1.0063 

3.1250 

-2.0969 

15 

.65 

1.0108 

3.2309 

-2.0933 

16 

.70 

1.0158 

3.3373 

-2.0896 

17 

.75 

1.0212 

3.4443 

-2.0856 

18 

.79 

1.0271 

3.5519 

-2.0814 

19 

.84 

1.0334 

3.6602 

-2,0770 

20 

.89 

1.0402 

3.7692 

-2.0723 


180 

8.34 

1.0632 

-.4841 

-2.0782 

181 

8.39 

1.0559 

-.3694 

-2.0828 

182 

8.44 

1.0490 

-.2557 

-2 . 087 1 

183 

8.48 

1.0424 

-.1430 

-2.0912 

184 

8.53 

1.0361 

-.0312 

-2.0950 

185 

8.58 

1.0302 

.0798 

-2.0986 

186 

8.62 

1.0247 

.1900 

-2.1019 

187 

8.67 

1.0195 

.2994 

-2.1049 

188 

8.72 

1 .0146 

.4080 

-2 . 1 077 

189 

8.76 

1.0101 

.5161 

-2,1103 

190 

8.81 

1.0059 

.6234 

-2.1126 

191 

8.86 

1.0020 

.7302 

-2.1146 

192 

8.90 

.9984 

.8365 

-2.1164 

193 

8.95 

.9952 

.9422 

-2.1180 

194 

9.00 

.9924 

1.0475 

-2.1193 

195 

9.04 

.9898 

1.1524 

-2.1204 

196 

9.09 

.9876 

1.2569 

-2.1212 

197 

9.14 

.9857 

1.3611 

-2.1218 

198 

9.18 

.9841 

1.4650 

-2.1222 

199 

9.23 

.9829 

1.5687 

-2,1223 

200 

9.28 

.9820 

1.6722 

-2. 1222 

201 

9.32 

.9814 

1.7755 

-2. 1219 

127C 

15.20 

2.0874 

-9.3223 

-1,3504 


COMPONENTS 


EXAMPLE 1 


)MBINED 

FAR FIELD 

DRAG NOISE — 
NEAR FIELD 

COMBINED 

OVERALL 

NOISE 

POIl 

NUMBl 

CPA) 

(PA) 

(PA) 

(PA) 

(PA) 


-.3464 

0.0000 

0.0000 

0.0000 

.6350 

1 

-.2427 

0.0000 

0.0000 

0.0000 

.7386 

2 

-.1387 

0.0000 

0.0000 

0.0000 

.8426 

3 

-.0346 

0.0000 

0.0000 

0.0000 

.9472 

4 

.0698 

0.0000 

0.0000 

0.0000 

1.0524 

5 

.1745 

0.0000 

0.0000 

0.0000 

1.1583 

6 

.2796 

0.0000 

0.0000 

0.0000 

1.2648 

7 

.3850 

0.0000 

0.0000 

0.0000 

1.3721 

8 

.4909 

0.0000 

0.0000 

0.0000 

1.4802 

9 

.5972 

0.0000 

0.0000 

0-0000 

1.5892 

10 

.7040 

0.0000 

0.0000 

0.0000 

1.6990 

11 

.8114 

0.0000 

0.0000 

0.0000 

1.8098 

12 

.9195 

0.0000 

0.0000 

0.0000 

1.9216 

13 

1.0281 

0.0000 

0.0000 

0.0000 

2.0344 

14 

1.1375 

0.0000 

0.0000 

0.0000 

2.1484 

15 

1.2477 

0.0000 

0.0000 

0.0000 

2.2635 

16 

1.3586 

0.0000 

0.0000 

0.0000 

2.3799 

17 

1.4705 

0.0000 

0.0000 

0.0000 

2.4976 

18 

1.5832 

0.0000 

0.0000 

0.0000 

2.6166 

19 

1.6969 

0.0000 

0.0000 

0.0000 

2.7370 

20 


-2.5622 

0.0000 

0.0000 

0.0000 

-1.4990 

180 


-2.4522 

0.0000 

0.0000 

0.0000 

-1.3963 

181 


-2.3428 

0.0000 

0.0000 

0.0000 

-1.2939 

182 


-2.2342 

0.0000 

0.0000 

0.0000 

-1.1918 

183 


-2.1262 

0.0000 

0.0000 

0.0000 

-1.0901 

184 


-2.0188 

0.0000 

0.0000 

0.0000 

-.9885 

185 


-1.9119 

0.0000 

0.0000 

0.0000 

-.8872 

186 


-1.8056 

0.0000 

0.0000 

0.0000 

-.7861 

187 


-1-6997 

0.0000 

0.0000 

0.0000 

-.6851 

188 


-1 .5942 

0.0000 

0.0000 

0.0000 

-.5842 

189 


-1.4891 

0.0000 

0.0000 

0.0000 

-.4833 

190 


-1 .3844 

0.0000 

0.0000 

0.0000 

-.3824 

191 


-1.2800 

0.0000 

0.0000 

0.0000 

-.2815 

192 


-1.1758 

0.0000 

0.0000 

0.0000 

-.1805 

193 


-1 .0718 

0.0000 

0.0000 

0.0000 

-.0794 

194 


-.9680 

0.0000 

0.0000 

0.0000 

.0218 

195 


-.8643 

0.0000 

0.0000 

0.0000 

.1233 

196 


-.7607 

0.0000 

0.0000 

0.0000 

.2250 

197 


-.6572 

0.0000 

0.0000 

0.0000 

.3269 

198 

W 

-.5537 

0.0000 

0.0000 

0.0000 

.4292 

199 

X 

-.4501 

0.0000 

0.0000 

0.0000 

.5319 

200 


-.3464 

0.0000 

0.0000 

0.0000 

.6350 

201 

M 

-10.6727 

0.0000 

0.0000 

0.0000 

-8.5853 
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APPENDIX C 


PRESSURE SPECTRA OF NOISE COMPONENTS 


EXAMPLE 1 


HARMONIC 

THICKNESS 

NUMBER 

FREQUENCY 

NOISE 


(HZ) 

(DB) 

1 

107.25 

99.01 

2 

214.50 

100.96 

3 

321.75 

100.96 

4 

429.00 

100.28 

5 

536.25 

99.26 

6 

643.50 

98.05 

7 

750.75 

96.73 

8 

858.00 

95.33 

9 

965.25 

93.89 

10 

1072.50 

92.42 

11 

1179.75 

90.91 

12 

1287.00 

89.37 

13 

1394.25 

87.83 

14 

1501 .50 

86.29 

15 

1608.75 

84.77 

16 

1716.00 

83.23 

17 

1823.25 

81.66 

18 

1930.50 

80.06 

19 

2037.75 

78.48 

20 

2145.00 

76.95 

21 

2252.25 

75.43 

22 

2359.50 

73.88 

23 

2466.75 

72.27 

24 

2574.00 

70.66 

25 

2681.25 

69.10 

26 

2788.50 

67.62 

27 

2895.75 

66.13 

28 

3003.00 

64.57 

29 

3110.25 

62.95 

30 

3217.50 

61.35 


LOADING 

DRAG 

OVERALL 

NOISE 

NOISE 

NOISE 

(DB) 

(DB) 

(DB) 

109.51 

-158.57 

109.88 

105.20 

-158.57 

106.69 

101.26 

-158.57 

104.24 

97.71 

-158.57 

102.30 

94.45 

-158.57 

100.58 

91.42 

-158.57 

98.96 

88.56 

-158.57 

97.38 

85.84 

-158.57 

95.81 

83.23 

-158.57 

94.25 

80.72 

-158.57 

92.68 

78.28 

-158.57 

91.11 

75.89 

-158.57 

89.52 

73.57 

-158.57 

87.94 

71.32 

-158.57 

86.37 

69.14 

-158.57 

84.82 

66.99 

-158.57 

83.26 

64.85 

-158.57 

81.67 

62.71 

-158.57 

80.06 

60.62 

-158.57 

78.47 

58.60 

-158.57 

76.92 

56.64 

-158.57 

75.40 

54.68 

-158.57 

73.84 

52.69 

-158.57 

72.23 

50.70 

-158.57 

70.60 

48.78 

-158.57 

69.04 

46.95 

-158.57 

67.55 

45.13 

-158.57 

66.06 

43.28 

-158.57 

64.50 

41 .41 

-158.57 

62.88 

39.63 

-158.57 

61.28 
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APPENDIX C 


SUBROUTINE FUNE2 

$ ( E2 , R , AA , AAP , CH , CHP , LED , LEDP , THKR AT , THKRATP , CAMRAT , CAMRATP ) 

REAL LED, LEDP 
ER = E2/R 

- TWIN OTTER INPUT FUNCTIONS, EXAMPLE 2. 

AAEXP = EXP(-3.491»ER) 

AA = 12.622 + 78.037*AAEXP 
AAP = -209.529*AAEXP 

ROOT r SQRT( 1.-ER ) 

POLY = 1.474 + ER»( 2.544 + ER ) 

CH = 0.069*ROOT»POLY 

CHP = 0.069/R*(-0.5/ROOT»POLY + ROOT»(2.»ER+2.544) ) 

LED = -0.25»CH 
LEDP = -0.25*CHP 

THKEXP = EXP(-11.200»ER) 

THKRAT = 0.069 + 3 .2244»THKEXP 
THKRATP = -27.778»THKEXP 

CAMRAT =0.03 
CAMRATP = 0. 

RETURN 

END 


SUBROUTINE FUNE2Q 

$ ( E2, R, Q, CMBR, CMBRP, THK, THKP ) 

- TWIN OTTER INPUT FUNCTIONS, EXAMPLE 2. 

CMBR = 4.»(Q-Q»»2) 

CMBRP = 4.-8.»Q 
C 

IF (Q.GT.0.5) GO TO 100 

THK = 0.989665*SQRT(Q) + Q»(-0. 23925 + Q»(-0.041 - Q»0.5594)) 
THKP = 0.4948325/SQRT(Q) - .23925 + Q»(-0.082 - Q»1.6782) 

GO TO 200 
C 

100 QM = 1 . - Q 
QMP = -1. 

THK = 0.01 + QM»(2.325 + QM»(-3.42 + QM*1.46)) 

THKP = (2.325 + QM»(-6.84 + QM»4.38) )»QMP 
C 

200 RETURN 
END 


EXAMPLE 2 
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APPENDIX C 


SUBROUTINE FUNPRES EXAMPLE 2 

$ (E2, R, Q, SURFP, DP, SPU, SPL, SIGMAU, SIGNAL) 

C 

LOGICAL SURFP 
C 

C - TWIN OTTER INPUT FUNCTIONS, EXAMPLE 2. 

C - SURFACE PRESSURE. 

C 

SURFP = .TRUE. 

A = . 7022^14*10. »»6. 

C 

ER = E2/R 

IF (Q-0.05) 10,20,20 
10 SPQU = EXP(800.»Q*(Q-.D) - 0.5 

SPQL = 0.9*EXP(500.»Q»(Q-.1)) - 0.4 
GO TO 100 
C 

20 IF ((Q-.05)»(Q-.5)) 30,30,40 

30 SPQL = -0.1421 

SPQU r -0.3647 
GO TO 100 
C 

40 SPQL = .2842»(Q-1 .0) 

SPQU = .7293*(Q-1 .0) 

C 

100 SPU = SPQU*ER»*2.833*SQRT(1 .-ER)»A 
SPL = SPQL»ER»»2.833*SQRT(1 .-ER)»A 
C 

DP = SPL - SPU 
C 

PI = 3.14159 
REV = 2145.0 
OMEGA = REV»PI/30. 

V3 = 63.79 

SIGMAU = 0.0018»((E2*OMEGA)»*2 + V3»V3) 

SIGNAL = SIGMAU 
C 

RETURN 

END 
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APPENDIX C 


DATA SHEET 


c 

= 

344.7 

M/SEC 




RHO 

= 

1.2060 KG/M»»3 




XO 

= 

( 7.28 

0.00, - 3 . 

.05) 

M OBSERVER IS IN MOTION 

V3 

r 

63.8 

M/SEC 





= 

142.7 

MPH 




R 

— 

1.300 

1 M 


HELICAL M-TIP = .867 

RINNER 

= 

.300 

M 




REV 

- 

2145.0 

RPM 


PER = 

9.3240E-03 SEC 

NBLADES 

= 

3 




BPF = 107.3 HZ 

NPCA 

:: 

15 




NPCAF = 3 

NLE 

= 

8 




NCHF = 2 

NTE 

= 

12 





THKMAXL 

= 

in 

o 

• 




NTAU = 10 

NPTS 

- 

200 



DT = 

4.6620E-05 SEC 

NSPEC 

= 

30 




FRACDT = .500 







MTRANS = .950 

EPSILON 

= 

.50 

% 



TRANS = .980 







INDEXSL = 0 

TORQUE 

r 

582.04 

N-M/BLADE 




THRUST 

= 

1894.25 

N/BLADE 




POWER 

= 

130.74 

KW/BLADE 






175.26 

HP/BLADE 




OASPL 

- 

115.83 

DB RE 20.E-6 

PA 





100.91 

DB RE 20.E-6 

PA 

(THICKNESS) 



= 

115.69 

DB RE 20.E-6 

PA 

(LOADING) 



= 

81 .64 

DB RE 20.E-6 

PA 

(DRAG) 



EXAMPLE 2 
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PRESSURE 

SIGNATURES OF 

NOISE COMPONENTS 



EXAMPLE 2 

POINT 

OBSERVER 

THICKNESS 


LOADING NOISE 






POINT 

NUMBER 








TIME 

NOISE 

FAR FIELD 

NEAR FIELD 

COMBINED 

FAR FIELD 

NEAR FIELD 

COMBINED 

NOISE 

NUMBER 


(MSEC) 

(PA) 

(PA) 

(PA) 

(PA) 

(PA) 

(PA) 

(PA) 

(PA) 

1 

0.00 

.7506 

-2.1545 

1.1437 

-1 .0108 

-.0402 

-.0195 

-.0596 

-.3199 

1 

2 

.05 

.7472 

-1.9772 

1.1337 

-.8435 

-.0370 

-.0196 

-.0566 

-.1530 

2 

3 

.09 

.7439 

-1 .8010 

1.1241 

-.6769 

-.0339 

-.0197 

-.0536 

.0134 

3 

4 

.14 

.7408 

-1.6258 

1.1149 

-.5109 

-.0308 

-.0198 

-.0506 

.1793 

4 

5 

.19 

.7379 

-1.4515 

1.1061 

-.3454 

-.0277 

-.0199 

-.0476 

.3449 

5 

6 

.23 

.7352 

-1 .2781 

1.0977 

-.1804 

-.0246 

-.0200 

-.0446 

.5102 

6 

7 

.28 

.7326 

-1 . 1054 

1 . 0897 

-.0157 

-.0216 

-.0201 

-.0416 

.6753 

7 

8 

.33 

.7303 

-.9335 

1.0821 

.1486 

-.0185 

-.0201 

-.0387 

.8402 

8 

9 

.37 

.7281 

-.7622 

1,0749 

.3128 

-.0155 

-.0202 

-.0358 

1.0051 

9 

10 

.42 

.7261 

-.5914 

1.0681 

.4767 

-.0125 

-.0203 

-.0328 

1.1700 

10 

11 

.47 

.7243 

-.4211 

1.0617 

.6406 

-.0096 

-.0203 

-.0299 

1.3350 

11 

12 

.51 

.7226 

-.2512 

1.0557 

.8045 

-.0066 

-.0204 

-.0270 

1.5001 

12 

13 

.56 

.7211 

-.0816 

1.0501 

.9685 

-.0037 

-.0204 

-.0241 

1.6656 

13 

14 

.61 

.7199 

.0878 

1.0449 

1.1326 

-.0007 

-.0205 

-.0212 

1.8313 

14 

15 

.65 

.7188 

.2570 

1.0400 

1.2970 

.0022 

-.0205 

-.0183 

1.9974 

15 

16 

.70 

.7179 

.4261 

1.0356 

1 .4616 

.0051 

-.0205 

-.0154 

2.1641 

16 

17 

.75 

.7171 

.5951 

1.0315 

1.6267 

.0080 

-.0205 

-.0125 

2.3313 

17 

18 

.79 

.7166 

.7643 

1.0279 

1.7922 

.0109 

-.0205 

-.0096 

2.4991 

18 

19 

.84 

.7163 

.9336 

1.0246 

1.9582 

.0138 

-.0205 

-.0067 

2.6677 

19 

20 

.89 

.7161 

1.1031 

1.0217 

2.1248 

.0167 

-.0205 

-.0038 

2.8371 

20 


180 

8.34 

181 

8.39 

182 

8.44 

183 

8.48 

184 

8.53 

185 

8.58 

186 

8.62 

187 

8.67 

188 

8.72 

189 

8.76 

190 

8.81 

191 

8.86 

192 

8.90 

193 

8.95 

194 

9.00 

195 

9.04 

196 

9.09 

197 

9.14 

198 

9.18 

199 

9.23 

200 

9.28 

201 

9.32 

139c 

15.76 


.8633 

-6.2696 

1.4537 

-4.8159 

-.1169 

-.0150 

-.1320 

-4.0846 

180 

,8562 

-6,0509 

1.4344 

-4.6165 

-.1127 

-.0153 

-.1280 

-3.8883 

181 

,8493 

-5.8350 

1.4155 

-4.4195 

-.1085 

-.0156 

-.1241 

-3.6943 

182 

8426 

-5.6218 

1.3971 

-4.2247 

-.1044 

-.0159 

-.1203 

-3.5024 

183 

.8360 

-5.4113 

1.3792 

-4.0321 

-.1004 

-.0162 

-.1166 

-3.3126 

184 

,8296 

-5.2033 

1.3618 

-3.8415 

-.0964 

-.0164 

-.1129 

-3.1248 

185 

,8233 

-4.9978 

1.3448 

-3.6530 

-.0925 

-.0167 

-.1092 

-2.9389 

186 

.8173 

-4.7947 

1.3283 

-3.4664 

-.0887 

-.0169 

-.1056 

-2.7547 

187 

,8114 

-4.5939 

1.3123 

-3.2816 

-.0849 

-.0172 

-.1021 

-2.5723 

188 

8056 

-4.3953 

1.2967 

-3.0986 

-.0812 

-.0174 

-.0986 

-2.3915 

189 

,8001 

-4.1988 

1.2816 

-2.9172 

-.0775 

-.0176 

-.0951 

-2.2123 

190 

,7947 

-4.0044 

1.2669 

-2.7375 

-.0739 

-.0178 

-.0917 

-2.0345 

191 

,7895 

-3.8119 

1.2526 

-2.5593 

-.0703 

-.0180 

-.0884 

-1 .8582 

192 

.7845 

-3.6213 

1.2388 

-2.3825 

-.0668 

-.0182 

-.0850 

-1.6831 

193 

,7796 

-3.4325 

1.2254 

-2 . 207 1 

-.0634 

-.0184 

-.0818 

-1.5092 

194 

,7749 

-3.2454 

1.2125 

-2.0330 

-.0599 

-.0186 

-.0785 

-1.3365 

195 

,7704 

-3.0600 

1 .2000 

-1 .8600 

-.0565 

-.0187 

-.0753 

-1.1649 

196 

,7661 

-2.8761 

1.1879 

-1.6882 

-.0532 

-.0189 

-.0721 

-.9942 

197 

,7620 

-2.6937 

1.1762 

-1.5175 

-.0499 

-.0191 

-.0689 

-.8245 

198 

,7580 

-2.5127 

1.1650 

-1.3477 

-.0466 

-.0192 

-.0658 

-.6556 

199 

7542 

-2.3330 

1.1541 

-1.1789 

-.0434 

-.0193 

-.0627 

-.4874 

200 

,7506 

-2.1545 

1.1437 

-1 .0108 

-.0402 

-.0195 

-.0596 

-.3199 

201 

6535 

-17.6574 

2.7168 

-14.9406 

-.3920 

.0082 

-.3837 

-14.6708 

139c 
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APPENDIX C 


PRESSURE SPECTRA OF NOISE COMPONENTS 


HARMONIC 

THICKNESS 

LOADING 

DRAG 

NUMBER 

FREQUENCY 

NOISE 

NOISE 

NOISE 


(HZ) 

(DB) 

(DB) 

(DB) 

1 

107.25 

95.00 

113.77 

79.24 

2 

214.50 

95.34 

109.27 

75.58 

3 

321.75 

93.83 

104.73 

71.70 

4 

429.00 

91.70 

100.38 

67.91 

5 

536.25 

89.28 

96.21 

64.23 

6 

643.50 

86.72 

92.20 

60.64 

7 

750.75 

84.07 

88.33 

57.14 

8 

858.00 

81.37 

84.54 

53.69 

9 

965.25 

78.59 

80.82 

50.28 

10 

1072.50 

75.79 

77.16 

46.91 

11 

1179.75 

73.01 

73.60 

43.63 

12 

1287.00 

70.29 

70.16 

40.42 

13 

1394.25 

67.54 

66.76 

37.20 

14 

1501.50 

64.69 

63.30 

33.92 

15 

1608.75 

61.78 

59.82 

30.66 

16 

1716.00 

59.02 

56.51 

27.55 

17 

1823.25 

56.49 

53.44 

24.60 

18 

1930.50 

53.95 

50.36 

21.57 

19 

2037.75 

51.06 

47.00 

18.29 

20 

2145.00 

47.82 

43.44 

14.88 

21 

2252.25 

44.75 

40.20 

11.76 

22 

2359.50 

42.30 

37.59 

9.15 

23 

2466.75 

40.33 

35.26 

6.72 

24 

2574.00 

38.20 

32.55 

4.01 

25 

2681.25 

35.63 

29.14 

.83 

26 

2788.50 

32.73 

25.33 

-2.62 

27 

2895.75 

29.53 

21.88 

-5.71 

28 

3003.00 

26.72 

20.44 

-7.36 

29 

3110.25 

25.97 

21.19 

-7.60 

30 

3217.50 

26.57 

21.13 

-8.63 


EXAMPLE 2 


OVERALL 

NOISE 

(DB) 

113.83 

109.44 

105.06 

100.93 

97.02 

93.30 

89.74 

86.29 
82.92 
79.62 
76.42 

73.35 

70.30 
67.20 
64.08 
61.12 
58.41 
55.71 

52.68 

49.36 

46.26 
43.78 

41.68 

39.27 

36.36 
33.39 
30.76 
29.05 
28.81 
28.85 
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SUBROUTINE FUNE2 

$ ( E 2 , R , AA , AAP , CH , CHP , LED , LEDP , THKRAT , THKRATP , CAMRAT , CAMRATP ) 

C 

REAL LED, LEDP, LED2 

DIMENSION ER1( 12) ,AA1( 12) ,ER2(20) ,CH2(20) ,LED2(20) ,THKRAT2(20) , 

$ W0RK1(60) ,WORK2(60) ,WORK3(60) ,WORK4(60) 

C 

C - HIGH SPEED PROPELLER INPUT FUNCTIONS, EXAMPLE 3. 

C 

C - DATA ARRAYS DEFINING BLADE GEOMETRY 
C AAl = BLADE ELEMENT ANGLE (DEG.) 

C CH2 = LOCAL CHORD DIVIDED BY PROPELLER DIAMETER 

LED2 = LEADING EDGE DISPLACEMENT DIVIDED BY LOCAL CHORD 
THKRAT2 = LOCAL THICKNESS RATIO 
C 

DATA ERl /. 239, .286, .367, .449, .531, .612, .694, .776, .857, 

$ .939, .980,1.000/ 

DATA AAl 782.83,80.65,76.99,73.61,70.21,66.89,63.75,60.80, 

$ 58.05,55.59,54.58,54.09/ 

DATA ER2/. 30,. 40, .50,. 55, .60, .65, .70, .75, .80, .825, .85, .875, .90, 

$ .925,. 95,. 96,. 97,. 98,. 99', 1.00/ 

DATA CH2 /. 170,. 187,. 200,. 202,. 200,. 195,. 188,. 178,. 166,. 160,. 152, 

$ .141, .129, .116, .102, .096, .090, .084, .077, .071 / 

DATA LED2 / . 527 , . 597 , .621 , . 598 , . 558 , .503 , .429, .332 , .206 , . 132 , .041 , 
$ -.075, -.218, -.399, -.637, -.760, -.896, -1.057, -1.260, -1.479/ 

DATA THKRAT2 /. 102, .060, .044, .039, .034, .029, .026, .024, .022, .022, 

$ . 021 , . 021 , . 021 , . 020 , . 020 , . 020 , . 020 , . 020 , . 020 , . 020 / 

DATA 11,12,13,14/-!, -1,-1,-!/ 

C 

ER = E2/R 
C 

OFFSET = -.435 

CALL PCS (ER1,AA1,12,ER,AA,AAP,I1,W0RK1) 

AA = AA + OFFSET 
AAP = AAP/R 
C 

IF (ER.LT.0.30) ER = 0.30 

CALL PCS (ER2,CH2,20,ER,CH,CHP,I2,W0RK2) 

CH = 2.*R*CH 
CHP = 2.*CHP 
C 

CALL PCS (ER2,LED2,20,ER,LED,LEDP,I3,W0RK3) 

LED = -CH*LED 

LEDP = -CHP*LED - CH*LEDP/R 
C 

CALL PCS (ER2,THKRAT2, 20, ER, THKRAT, THKRATP, 14, WORK4) 

THKRATP = THKRATP/ R 
C 

CAMRAT = 0.03 - 0.0913*(E2 - 0.093) 

CAMRATP = -0.0913 
C 

RETURN 

END 
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SUBROUTINE FUNE2Q 

$ ( E2, R, Q, CMBR, CMBRP, THK, THKP ) 

- HIGH SPEED PROPELLER INPUT FUNCTIONS, EXAMPLE 3. 

ER = E2/R 

CMBR = 4.»(Q-Q»»2) 

CMBRP = 4.-8.»Q 

IF (ER.GT.0.45) GO TO 50 

- SERIES 65 AIRFOIL. 

IF (Q.GT.0.4) GO TO 25 

THK = 1.1722»SQRT(Q) + Q«(0.078 + Q»(-12.3902 + 
$ Q*(56.9872 - 0*75.6761))) 

THKP = 0.5861/SQRT(Q) + 0.078 + Q*(-24.7804 + 

$ Q*(170.9616 - 0*302.7044)) 

GO TO 200 


OM = 0 - 0.4 

THK = 0.5 + 0M*0M*(-2.8333 + 0M*2.4074) 

THKP = 0M*(-5.6666 + 0M*7.2222) 

GO TO 200 

- SERIES 16 AIRFOIL. 

50 IF (O.GT.0.5) GO TO 100 

THK = 0.989665*S0RT(0) + 0*(-0. 23925 + 0*(-0.041 - 0*0.5594)) 
THKP = 0.4948325/S0RT(0) - .23925 + 0*(-0.082 - 0*1.6782) 

GO TO 200 
C 

100 OM = 1. - 0 
OMP = -1. 

THK = 0.01 + 0M*(2.325 + 0M*(-3.42 + 0M*1.46)) 

THKP = (2.325 + 0M*(-6.84 + OM*4.38))*OMP 
C 

200 RETURN 
C 

END 
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SUBROUTINE FUNPRES 

$ ( E2 , R , Q , SURFP , DP , SPU , SPL , SIGMAU , SIGNAL ) 

LOGICAL SURFP 

DIMENSION ER2(20), CP2(20), W0RK(60) 

C 

C - HIGH SPEED PROPELLER INPUT FUNCTIONS, EXAMPLE 3. 

C - DIFFERENTIAL PRESSURE INPUT. 

C 

DATA ER2/. 30,. 40,. 50,. 55,. 60,. 65,. 70,. 75,. 80,. 825,. 85,. 875,. 90, 
$ .925,.95,.96,.97,.98,.99,1.00/ 

DATA CP2 /. 214, .259, .307, .315, .318, .319, .320, .327, .339, .340, 

$ .336, .329, .320, .312, .297, .285, .278, .258, .248,0./ ' 

DATA INTRVL/-1/ 

C 

ER = E2/R 

IF (ER.LT.0.30) ER = 0.30 
SURFP = .FALSE. 

- DEFINE POWER FACTOR (FACTOR), ROTATIONAL SPEED (OMEGA), 

FORWARD SPEED (V3), DENSITY (RHO). 

FACTOR = 0.7861 
PI = 3.14159 
OMEGA = 7878. 4 / 30 . »PI 
V3 = 241.2 
RHO = .4570 

CALL PCS (ER2,CP2,20,ER,CP,CPP,INTRVL,W0RK) 

CP = FACTOR»CP 

- DEFINE LOCAL VELOCITY, SQUARED (VINES), AND 
DIFFERENTIAL PRESSURE (DP). 

VINFS = (E2*0MEGA)»»2. + V3*V3 
DP = 3.*RHO»VINFS»CP*(Q-Q»Q) 

C 

SPU = 0. 

SPL = 0. 

SIGMAU = 0. 

SIGNAL = 0. 

C 

RETURN 

END 
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DATA SHEET 


c 

= 

301.5 

M/SEC 



RHO 

= 

.4570 KG/M**3 



XO 

= 

( .81 

, 0.00, 

.01) M OBSERVER IS IN MOTION 

V3 

= 

241.2 

M/SEC 




= 

539.5 

MPH 



R 


.312 

M 

HELICAL M-TIP = 1.170 

RINNER 

= 

.093 

I M 



REV 

= 

7878.4 

RPM 

PER = 

1.9039E-03 SEC 

NBLADES 

= 

4 



BPF = 525.2 HZ 

NPCA 

= 

20 



NPCAF = 3 

NLE 

= 

5 



NCHF = 2 

NTE 

= 

10 




THKMAXL 

= 

.05 



NT AH =10 

NPTS 

= 

150 


DT = 

1.2693E-05 SEC 

NSPEC 

= 

30 



FRACDT = 2.000 






MTRANS = .950 

EPSILON 

= 

.50 

% 


TRANS = .980 






INDEXSL = 38794 

TORQUE 

= 

20.22 

N-M/ BLADE 



THRUST 

= 

46.51 

N/BLADE 



POWER 

= 

16.69 

KW/ BLADE 




= 

22.37 

HP/ BLADE 



OASPL 

= 

142.38 

DB RE 20.E-6 

PA 




139.78 

DB RE 20.E-6 

PA (THICKNESS) 



= 

133.45 

DB RE 20.E-6 

PA (LOADING) 



= 

-156.07 

DB RE 20.E-6 

PA (DRAG) 
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PRESSURE SIGNATURES OF NOISE COMPONENTS 


EXAMPLE 3 


POINT 

OBSERVER 

THICKNESS 


NUMBER 

TIME 

NOISE 

FAR FIELD 


(MSEC) 

(PA) 

(PA) 

1 

0.00 

-68.6105 

95.7354 

2 

.01 

-56.0640 

91.1677 

3 

.03 

-115.4197 

81.9249 

4 

.04 

-104.9222 

82.5168 

5 

.05 

-74.7874 

79.5310 

6 

.06 

-103.5242 

75.1177 

7 

.08 

-110.7724 

65.4773 

8 

.09 

-129.4907 

59.7124 

9 

.10 

-133.2506 

51.7046 

10 

.11 

-137.1576 

55.2064 

11 

.13 

-150.5800 

29.7830 

12 

.14 

-168.6043 

29.0084 

13 

.15 

-161.5778 

23.7662 

14 

.17 

-143.3184 

20.6980 

15 

.18 

-160.5617 

20.3922 

16 

.19 

-161.6439 

10.1030 

17 

.20 

-168.6064 

-11.8672 

18 

.22 

-179.0823 

-22.0207 

19 

.23 

-156.1153 

-24.6496 

20 

.24 

-169.9904 

-14.1905 


LOADING NOISE - 
NEAR FIELD 

COMBINED 

FAR FIELD 

(PA) 

(PA) 

(PA) 

-15.7404 

79.9950 

0.0000 

-15.5995 

75.5681 

0.0000 

-15.5212 

66.4037 

0.0000 

-15.4019 

67.1149 

0.0000 

-15.3685 

64.1625 

0.0000 

-15.2701 

59.8476 

0.0000 

-15.2323 

50.2450 

0.0000 

-15.1169 

44.5956 

0.0000 

-15.1019 

36.6026 

0.0000 

-15.0537 

40.1527 

0.0000 

-15.0287 

14.7543 

0.0000 

-15.0054 

14.0029 

0.0000 

-14.9731 

8.7931 

0.0000 

-14.9797 

5.7183 

0.0000 

-14.9839 

5.4083 

0.0000 

-15.0199 

-4.9169 

0.0000 

-14.9720 

-26.8392 

0.0000 

-15.0324 

-37.0531 

0.0000 

-15.0364 

-39.6860 

0.0000 

-15.2049 

-29.3954 

0.0000 


DRAG NOISE OVERALL POINT 

NEAR FIELD COMBINED NOISE NUMBER 


(PA) 

(PA) 

(PA) 


0.0000 

0.0000 

11.3845 

1 

0.0000 

0.0000 

19.5041 

2 

0.0000 

0.0000 

-49.0160 

3 

0.0000 

0.0000 

-37.8073 

4 

0.0000 

0.0000 

-10.6249 

5 

0.0000 

0.0000 

-43.6766 

6 

0.0000 

0.0000 

-60.5274 

7 

0.0000 

0.0000 

-84.8951 

8 

0.0000 

0.0000 

-96.6479 

9 

0.0000 

0.0000 

-97.0048 

10 

0.0000 

0.0000 

-135.8257 

11 

0.0000 

0.0000 

-154.6014 

12 

0.0000 

0.0000 

-152.7847 

13 

0.0000 

0.0000 

-137.6001 

14 

0.0000 

0.0000 

-155.1534 

15 

0.0000 

0.0000 

-166.5608 

16 

0.0000 

0.0000 

-195.4456 

17 

0.0000 

0.0000 

-216.1354 

18 

0.0000 

0.0000 

-195.8013 

19 

0.0000 

0.0000 

-199.3858 

20 


O 


1^0 1.76 -46.9409 129.7920 -17.2485 112.5435 0.0000 

1^1 1.78 -24.3012 134.1412 -17.0903 117.0509 0.0000 

1^2 1.79 -40.0762 127.8688 -16.9421 110.9267 0.0000 

1^3 1.80 -41.5873 127.3471 -16.7731 110.5740 0.0000 

1^^ 1.82 -65.2115 118.5694 -16.6325 101.9369 0.0000 

1^5 1.83 -61.3185 120.5018 -16.4795 104.0223 0.0000 

1^6 1.84 -67.9980 113.1377 -16.3405 96.7972 0.0000 

147 1.85 -69.5087 115.5055 -16.2138 99.2917 0.0000 

148 1.87 -61.9047 103.0294 -16.0586 86.9709 0.0000 

149 1.88 -55.9023 107.9377 -15.9703 91.9674 0.0000 

150 1.89 -58.6788 102.5642 -15.8180 86.7462 0.0000 

151 1.90 -68.6105 95.7354 -15.7404 79.9950 0.0000 

109C 3.27 183.7243 75.7185 -21.5544 54.1641 0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


65.6025 140 

92.7497 141 

70.8504 142 

68.9867 143 

36.7254 144 

42.7038 145 

28.7992 146 

29.7830 147 

25.0661 148 

36.0651 149 

28.0674 150 

11.3845 151 

M 

237.8883 109C ^ 
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EXAMPLE 3 

PRESSURE SPECTRA OF NOISE COMPONENTS 


— 

HARMONIC 

THICKNESS 

LOADING 

DRAG 

OVERALL 

NUMBER FREQUENCY 

NOISE 

NOISE 

NOISE 

NOISE 


(HZ) 

(DB) 

(DB) 

(DB) 

(DB) 

1 

525.23 

138.48 

132.83 

-156.07 

141.34 

2 

1050.45 

128.37 

124.31 

-156.07 

132.30 

3 

1575.68 

127.98 

112.38 

-156.07 

129.24 

4 

2100.91 

124.73 

105.41 

-156.07 

124.46 

5 

2626.13 

119.54 

101.35 

-156.07 

119.99 

6 

3151.36 

123.16 

97.07 

-156.07 

123.54 

7 

3676.59 

119.60 

96.02 

-156.07 

119.85 

8 

4201.81 

121.17 

88.56 

-156.07 

121.25 

9 

4727.04 

119.50 

95.93 

-156.07 

119.42 

10 

5252.27 

117.58 

96.89 

-156.07 

116.96 

11 

5777.49 

118.10 

90.79 

-156.07 

117.93 

12 

6302.72 

115.86 

93.44 

-156.07 

115.86 

13 

6827.95 

115.04 

88.17 

-156.07 

114.95 

14 

7353.17 

111.45 

90.78 

-156.07 

111.27 

15 

7878.40 

108.66 

88.41 

-156.07 

108.23 

16 

8403.63 

110.01 

89.78 

-156.07 

109.36 

17 

8928.85 

95.37 

91.75 

-156.07 

86.72 

18 

9454.08 

101.08 

77.83 

-156.07 

100.86 

19 

9979.31 

95.67 

90.49 

-156.07 

99.12 

20 

10504.53 

93.05 

74.55 

-156.07 

93.23 

21 

11029.76 

90.89 

89.30 

-156.07 

89.44 

22 

11554.99 

96.16 

84.46 

-156.07 

98.12 

23 

12080.21 

84.11 

88.56 

-156.07 

91.98 

24 

12605.44 

101.14 

90.26 

-156.07 

101.53 

25 

13130.67 

70.00 

85.25 

-156.07 

83.67 

26 

13655.89 

91.63 

88.00 

-156.07 

89.16 

27 

14181.12 

98.15 

89.59 

-156.07 

100.70 

28 

14706.35 

84.62 

84.44 

-156.07 

87.66 

29 

15231.57 

103.21 

90.36 

-156.07 

104.92 

30 

15756.80 

102.52 

89.51 

-156.07 

102.03 
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SUBROUTINE FUNE2 

$ (E2 , R, AA , AAP , CH , CHP , LED , LEDP , THKRAT , THKRATP , CAMRAT , CAMRATP ) 

C 

REAL LED, LEDP, LED2 

DIMENSION ER1( 12) ,AA1(12) ,ER2(20) ,CH2(20) ,LED2(20) ,THKRAT2(20) , 

$ W0RKK60) ,WORK2(60) ,WORK3(60) ,WORK4(60) 

C 

C - HIGH SPEED PROPELLER INPUT FUNCTIONS, EXAMPLE 4. 

C 

C - DATA ARRAYS DEFINING BLADE GEOMETRY 
C AAl = BLADE ELEMENT ANGLE (DEG.) 

C CH2 = LOCAL CHORD DIVIDED BY PROPELLER DIAMETER 

C LED2 = LEADING EDGE DISPLACEMENT DIVIDED BY LOCAL CHORD 
C THKRAT2 = LOCAL THICKNESS RATIO 
C 

DATA ERl /. 239, .286, .367, .449, .531, .612, .694, .776, .857, 

$ .939, .980,1.000/ 

DATA AAl 782.83,80.65,76.99,73.61,70.21,66.89,63.75,60.80, 

$ 58.05,55.59,54.58,54.09/ 

DATA ER2/.30, .40, .50, .55, .60, .65, .70, . 75 , .80 , .825 , .85 , .875 ,. 90 , 

$ .925, .95, .96, .97, .98, .99, 1.00/ 

DATA CH2 /. 170,. 187,. 200,. 202,. 200,. 195,. 188,. 178,. 166,. 160,. 152, 

$ .141, .129, .116, .102, .096, .090, .084, .077, .071 / 

DATA LED2 /.527 , . 597 , .621 , .598 , . 558 , .503 , .429 , .332 , . 206 , . 132 , .041 , 
$ -.075, -.218, -.399, -.637, -.760, -.896, -1.057, -1.260, -1.479/ 

DATA THKRAT2 / . 102 , .060 , .044 , .039 , .034 , .029 , .026 , .024 , .022 , .022 , 

$ .021 , .021 , .021 , .020, .020, .020, .020, .020, .020 , .020/ 

DATA 11,12,13,14/-!, -1,-1,-!/ 

C 

ER = E2/R 
C 

OFFSET = -.435 

CALL PCS (ER1,AA1,12,ER,AA,AAP,I1,W0RK1) 

AA = AA + OFFSET 
AAP = AAP/R 
C 

IF (ER.LT.0.30) ER = 0.30 

CALL PCS (ER2,CH2,20,ER,CH,CHP,I2,WORK2) 

CH = 2.*R*CH 
CHP = 2.*CHP 
C 

CALL PCS (ER2,LED2,20,ER,LED,LEDP,I3,WORK3) 

LED = -CH*LED 

LEDP = -CHP*LED - CH*LEDP/R 
C 

CALL PCS (ER2,THKRAT2, 20, ER, THKRAT, THKRATP, 14, W0RK4) 

THKRATP = THKRATP /R 
C 

CAMRAT = 0.03 - 0.0913*(E2 - 0.093) 

CAMRATP = -0.0913 
C 

RETURN 

END 
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EXAMPLE 4 


SUBROUTINE FUNE2Q 

$ ( E2, R, Q, CMBR, CMBRP, THK, THKP ) 

- HIGH SPEED PROPELLER INPUT FUNCTIONS, EXAMPLE 4. 

ER = E2/R 

CMBR = 4.»(Q-Q»*2) 

CMBRP = 4.-8.*Q 

IF (ER.GT.0.45) GO TO 50 

- SERIES 65 AIRFOIL. 

IF (Q.GT.0.4) GO TO 25 

THK = 1 .1722«SQRT(Q) + Q*(0.078 + Q*(-12.3902 + 

$ Q*(56.9872 - 0*75.6761))) 

THKP = 0.5861/SQRT(Q) + 0.078 + Q*(-24.7804 + 

$ Q*(170.9616 - 0*302.7044)) 

GO TO 200 

; OM = 0 - 0.4 

THK = 0.5 + 0M*0M*(-2.8333 + 0M*2.4O74) 

THKP = 0M*(-5.6666 + 0M*7.2222) 

GO TO 200 

- SERIES 16 AIRFOIL. 

50 IF (O.GT.0.5) GO TO 100 

THK = 0.989665*S0RT(0) + 0*(-0. 23925 + 0*(-0.041 - 0*0.5594)) 
THKP = 0.4948325/S0RT(0) - .23925 + 0*(-0.082 - 0*1.6782) 

GO TO 200 
C 

100 OM = 1 . - 0 
OMP = -1. 

THK = 0.01 + 0M*(2.325 + 0M*(-3.42 + 0M*1.46)) 

THKP = (2.325 + 0M*(-6.84 + 0M*4 .38) ) *0MP 
C 

200 RETURN 
C 

END 
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SUBROUTINE FUNPRES 

$ ( E2 , R , Q , SURFP , DP , SPU , SPL , SIGMAU , SIGMAL) 

LOGICAL SURFP 

DIMENSION ER2(20), CP2(20), W0RK(60) 

- HIGH SPEED PROPELLER INPUT FUNCTIONS, EXAMPLE 4. 

- DIFFERENTIAL PRESSURE INPUT. 

DATA ER2/. 30,. 40,. 50,. 55,. 60,. 65,. 70,. 75,. 80,. 825,. 85,. 875,. 90, 
$ .925, .95, .96, .97, .98, .99,1.00/ 

DATA CP2 /. 214, .259, .307, .315, .318, .319, .320, .327, .339, .340, 

$ .336, .329, .320, .312, .297, .285, .278, .258, .248,0./ 

DATA INTRVL/-1/ 

ER = E2/R 

IF (ER.LT.0.30) ER = 0.30 
SURFP = .FALSE. 

- DEFINE POWER FACTOR (FACTOR) , ROTATIONAL SPEED (OMEGA) , 

FORWARD SPEED (V3) , DENSITY (RHO). 

FACTOR r 0.7861 
PI = 3.14159 
OMEGA = 7878. 4/30. «PI 
V3 = 241.2 
RHO = .4570 

CALL PCS (ER2,CP2,20,ER,CP,CPP,INTRVL,W0RK) 

CP = FACTOR»CP 

- DEFINE LOCAL VELOCITY, SQUARED (VINES), AND 
DIFFERENTIAL PRESSURE (DP). 

VINFS = (E2»0MEGA)*»2. + V3*V3 
DP = 3.*RHO»VINFS»CP*(Q-Q»Q) 

C 

SPU = 0. 

SPL = 0. 

SIGMAU = 0. 

SIGMAL = 0. 

C 

RETURN 

END 
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DATA SHEET 


C = 301.5 M/SEC 

RHO = .4570 KG/M**3 


XO = ( .81, 0.00, -.81) M OBSERVER IS IN MOTION 


V3 

= 

241.2 

M/SEC 





= 

539.5 

MPH 




R 

= 

.312 

, M 


HELICAL M-TIP = 1.170 

RINNER 

= 

.093 

M 




REV 

= 

7878.4 

RPM 


PER = 

1.9039E-03 SEC 

NBLADES 

= 

4 




BPF = 525.2 HZ 

NPCA 


20 




NPCAF = 1 

NLE 

= 

5 




NCHF = 1 

NTE 

= 

10 





THKMAXL 

= 

.05 




NTAU = 10 

NPTS 

= 

150 



DT = 

1.2693E-05 SEC 

NSPEC 

= 

30 




FRACDT = 1.000 







MTRANS = .950 

EPSILON 

= 

.50 

% 



TRANS = .980 







INDEXSL = 0 

TORQUE 

= 

20.36 

N-M/ BLADE 




THRUST 

= 

46.82 

N/ BLADE 




POWER 

= 

16.79 

KW/ BLADE 





= 

22.51 

HP /BLADE 




OASPL 

= 

124.48 

DB RE 20.E-6 

PA 




= 

104.57 

DB RE 20.E-6 

PA (THICKNESS) 



= 

125.17 

DB RE 20.E-6 

PA 

(LOADING) 



= 

-156.07 

DB RE 20.E-6 

PA 

(DRAG) 



EXAMPLE 4 
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U1 

to 


PRESSURE SIGNATURES OF NOISE COMPONENTS 


EXAMPLE 4 


POINT 

OBSERVER 

THICKNESS 


LOADING NOISE 

NUMBER 

TIME 

NOISE 

FAR FIELD 

NEAR FIELD 


(MSEC) 

(PA) 

(PA) 

(PA) 

1 

0.00 

4.1082 

-4.4962 

1.1994 

2 

.01 

4,0377 

-3.2728 

1.0313 

3 

.03 

3.9685 

-2.0617 

.8685 

4 

.04 

3.9013 

-.8638 

,7111 

5 

.05 

3.8345 

.3234 

.5590 

6 

.06 

3.7695 

1.4988 

.4125 

7 

.08 

3.7054 

2.6664 

.2713 

8 

.09 

3.6422 

3.8225 

.1356 

9 

.10 

3.5804 

4.9686 

.0055 

10 

.11 

3.5191 

6.1049 

-.1191 

11 

.13 

3.4590 

7.2328 

-.2382 

12 

.14 

3.3996 

8.3510 : 

-.3517 

13 

.15 

3.3410 

9.4610 

-.4595 

14 

.17 

3.2832 

10.5636 

-.5618 

15 

.18 

3,2259 

11.6572 

-.6583 

16 

.19 

3.1693 

12,7447 

-.7493 

17 

.20 

3.1132 

13.8240 

-.8345 

18 

.22 

3.0581 

14.8964 

-.9139 

19 

,23 

3.0028 

15.9630 

-.9877 

20 

.24 

2,9484 

17.0218 

-1.0557 


DRAG NOISE 


COMBINED 

FAR FIELD 

NEAR FIELD 

COMBINED 

(PA) 

(PA) 

(PA) 

(PA) 

-3.2967 

0.0000 

0.0000 

0.0000 

-2.2415 

0,0000 

0.0000 

0.0000 

-1.1932 

0.0000 

0.0000 

0.0000 

-.1527 

0.0000 

0.0000 

0.0000 

.8825 

0.0000 

0.0000 

0.0000 

1.9112 

0.0000 

0.0000 

0.0000 

2.9377 

0.0000 

0.0000 

0.0000 

3.9581 

0.0000 

0.0000 

0.0000 

4.9740 

0.0000 

0.0000 

0.0000 

5.9858 

0.0000 

0.0000 

0.0000 

6.9947 

0.0000 

0.0000 

0.0000 

7.9994 

0.0000 

0,0000 

0.0000 

9.0015 

0.0000 

0.0000 

0.0000 

10.0018 

0.0000 

0.0000 

0.0000 

10.9988 

0.0000 

0.0000 

0.0000 

11.9954 

0.0000 

0.0000 

0,0000 

12.9895 

0.0000 

0.0000 

0.0000 

13.9824 

0.0000 

0.0000 

0.0000 

14.9753 

0.0000 

0.0000 

0.0000 

15.9662 

0.0000 

0.0000 

0.0000 


OVERALL POINT 

NOISE NUMBER 

(PA) 


.8115 

1 

1.7963 

2 

2.7754 

3 

3.7486 

4 

4.7170 

5 

5.6807 

6 

6.6431 

7 

7.6003 

8 

8.5544 

9 

9.5048 

10 

10.4537 

11 

11.3990 

12 

12.3425 

13 

13.2850 

14 

14.2247 

15 

15.1648 

16 

16.1027 

17 

17.0405 

18 

17.9780 

19 

18.9146 

20 


140 

1.76 

141 

1.78 

142 

1.79 

143 

1.80 

144 

1.82 

145 

1.83 

146 

1.84 

147 

1.85 

148 

1.87 

149 

1.88 

150 

1.89 

151 

1.90 


5.0039 -18.8979 
4.9111 -17.5075 
4.8208 -16.1370 
4.7332 -14.7802 
4.6490 -13.4435 
4.5658 -12.1196 
4.4851 -10.8140 
4.4060 -9.5209 
4.3289 -8.2447 
4.2535 -6.9806 
4.1800 -5.7325 
4.1082 -4.4962 


6C 1.97 


3.7695 1.5012 


3.3900 

-15.5079 

0.0000 

0.0000 

0.0000 

-10.5040 

140 

3.1659 

-14.3415 

0.0000 

0.0000 

0.0000 

-9.4304 

141 

2.9467 

-13.1903 

0.0000 

0.0000 

0.0000 

-8.3696 

* 142 

2.7323 

-12.0480 

0,0000 

0.0000 

0.0000 

-7.3148 

143 

2.5228 

-10.9206 

0.0000 

0.0000 

0.0000 

-6.2716 

144 

2.3184 

-9.8012 

0.0000 

0.0000 

0.0000 

-5.2355 

145 

2.1190 

-8.6951 

0.0000 

0.0000 

0.0000 

-4.2100 

146 

1.9247 

-7.5962 

0.0000 

0.0000 

0.0000 

-3. 1902 

147 

1.7356 

-6.5092 

0.0000 

0.0000 

0.0000 

-2.1803 

148 

1.5516 

-5.4290 

0.0000 

0.0000 

0.0000 

-1.1755 

149 

1.3729 

-4.3596 

0.0000 

0.0000 

0.0000 

-.1796 

150 

1.1994 

-3.2967 

0,0000 

0.0000 

0.0000 

.8115 

151 

.4124 

1.9136 

0.0000 

0.0000 

0,0000 

5.6832 

6C 
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APPENDIX C 


PRESSURE SPECTRA OF NOISE COMPONENTS 


— 

HARMONIC 

THICKNESS 

LOADING 

DRAG 

OVERALL 

NUMBER FREQUENCY 

(HZ) 

NOISE 

(DB) 

NOISE 

(DB) 

NOISE 

(DB) 

NOISE 

(DB) 

1 

525.23 

102.13 

124.28 

-156.07 

123.64 

2 

1050.45 

98.55 

116.73 

-156.07 

115.83 

3 

1575.68 

94.74 

110.19 

-156.07 

109.05 

4 

2100.91 

90.98 

104.29 

-156.07 

102.92 

5 

2626.13 

87.24 

98.79 

-156.07 

97.20 

6 

3151.36 

83.48 

93.55 

-156.07 

91 .76 

7 

3676.59 

79.70 

88.48 

-156.07 

86.50 

8 

4201.81 

75.87 

83.52 

-156.07 

81 .38 

9 

4727.04 

71 .99 

78.64 

-156.07 

76.36 

10 

5252.27 

68.06 

73.80 

-156.07 

71 .42 

11 

5777.49 

64.09 

69.00 

-156.07 

66.57 

12 

6302.72 

60.06 

64.24 

-156.07 

61 .81 

13 

6827.95 

55.98 

59.47 

-156.07 

57.07 

14 

7353.17 

51 .86 

54.78 

-156.07 

52.52 

15 

7878.40 

47.81 

50.12 

-156.07 

47 ,. 85 

16 

8403.63 

43.69 

45.40 

-156.07 

43.40 

17 

8928.85 

39.55 

40.88 

-156.07 

39.10 

18 

9454.08 

35.50 

36.22 

-156.07 

34.24 

19 

9979.31 

31 .41 

30.40 

-156.07 

30.42 

20 

10504.53 

27.25 

26.33 

-156.07 

26.35 

21 

11029.76 

23.30 

20.95 

-156.07 

22.35 

22 

11554.99 

20.82 

15.54 

-156.07 

19.98 

23 

12080.21 

16.24 

18.56 

-156.07 

20.62 

24 

12605.44 

13.46 

14.94 

-156.07 

8.78 

25 

13130.67 

9.04 

3.50 

-156.07 

11.89 

26 

13655.89 

10.99 

-.02 

-156.07 

9.47 

27 

14181.12 

8.08 

8.62 

-156.07 

13.23 

28 

14706.35 

8.96 

- 2.83 

-156.07 

9.55 

29 

15231.57 

4.35 

14.08 

-156.07 

11.84 

30 

15756.80 

9.71 

16.95 

-156.07 

12.33 
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APPENDIX D 


EXAMPLES OF LANGLEY GRAPHICS OUTPUT 

The following graphs were selected from the graphics output generated by 
program NPROPFN for the four examples given in appendix C. Each graph is 
clearly labeled with an example number at the top right corner. The graphs 
which have been included for each example are: 

Overall Pressure 

Overall Spectrum 

Thickness Noise 

Thickness Spectrum 

Loading Noise Components 

Loading Noise 

Loading Spectrum 

For example 2 , three additional graphs have been included: 

Drag Noise Components 
Drag Noise 
Drag Spectrum 

The Langley version of this program also has the capability to plot the 
r\i and ri3 blade coordinates, angle of attack, chord, and thickness ratio 
functions as each varies with the blade radius 1)2 • 
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